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FOREWORD 


This volume presents the results of an analytical study 
performed from May 18, 1967 through May 17, 1968 on the use of 
dynamic scale models to determine launch vehicle characteristics. 
The investigation was conducted by the Denver Division of the 
Martin Marietta Corporation, Denver, Colorado, for the National 
Aeronautics and Space Administration, George C, Marshall Space 
flight Center, Huntsville, Alabama, under Contract No. NAS8-21101 . 
Mr, L. Kief ling was the principal representative for the contract- 
ing office. 

Mr. George Morosow was the Program Manager for the Denver 
Division and all work was performed under the direction of 
Mr. Morosow and Mr. Ivan J. Jaszlics, Principal Investigator, 
Significant contributions were provided by C, C. Feng, R. C, 
Reuter, Jr. and J. N. Singh. The report was critically reviewed 
by Professor L. G. Tulin of the University of Colorado, 

Results of an experimental study to confirm some of the 
analytical results presented herein are provided in Volume II, 
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ABSTRACT 


The purpose of this study was to develop modeling techniques 
that can be used for launch vehicle and space structure dynamic 
simulation when the parameters which define similarity between 
prototype and model cannot, due to physical or manufacturing 
limitations, be fully satisfied. In practice, all but the largest 
models have to make compromises in structural simulation and even 
very large models, using direct replica scaling, cannot duplicate 
all essential parameters. 

During the study it became obvious that there is need for a 
comprehensive summary of the scaling techniques applicable to the 
design of launch vehicle and space structure models. This report, 
provides that summary; describes general model scaling methods; 
developes basic similarity laws for various structures including 
liquid propellants and shell structures. A survey of material and 
manufacturing limitations and a summary of the economic aspects 
of model construction and testing is included. 
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SUMMARY 


Results of an analytical study on the use of dynamic scale models 
to determine launch vehicle characteristics are presented . The report 
is written in a way that permits its use as a general guide for 
launch vehicle model scaling and design. 

The dimensional analysis methods required to establish dynamic 
model scaling parameters are developed in matrix form and used to 
summarize the scaling laws applicable to launch vehicle structures. 
Scaling laws are presented for beam-like and shell-like structures. 
Visco-elastic effects and coupled structural-liquid slosh effects 
are considered. The scaling laws applicable to the simulation of 
liquid propellant dynamic phenomena are presented. Some special 
distorted modeling techniques for shells are developed and the 
necessity of an experimental program to confirm these techniques is 
emphasized. 

A limited discussion of the scaling requirements for space 
structures is presented and the relations for thermal scaling and 
scaling of gravity-gradient induced forces are summarized. A test 
configuration to simulate the gravity-gradient is discussed. A 
detailed theoretical and experimental program to investigate these 
and other modeling techniques applicable to large space structures 
is recommended. 

A summary of the economic aspects of dynamic model design and 
construction and a review of past and present launch vehicle dynamic 
models are included. 
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1. INTRODUCTION 


The first decade of space exploration was characterized by the 
continual increase in launch vehicle and payload dimensions and, with 
this increase, structural frequencies were reduced. Fundamental 
bending frequencies of first generation vehicles were on the order 
of 10 Hz; Titan and Saturn class vehicles have fundamental bending 
frequencies, during the boost phase, of 1 to 2 Hz. The degree of 
dynamic interaction between structural vibrations, propellant slosh 
and vehicle control systems increased considerably. 

All present generation launch vehicle structures are designed, 
at least in part, by the dynamic loading conditions that arise during 
the pre-launch or during the boost phase. These conditions include 
loads due to ground winds, lift-off loads, dynamic loads generated 
by wind shear or atmospheric turbulence, transonic buffet loads, 
staging loads, upper stage engine ignition loads and engine shutdown 
transient loads. 

The complete evaluation of these loads is only possible through 
combined analytical and test efforts. While analysis can predict the 
fundamental frequency of most conventional launch vehicle structures 
with approximately 3% to 8% accuracy, the risk involved in using only 
analytical methods, especially if the structure is complex, is so great 
that experimental verification is required. The analytical methods 
available for prediction of transonic buffet loads and ground wind- 
induced loads are of questionable accuracy at the present time: 
the use of test results in conjunction with analytical studies is of 
the greatest importance in predicting these aeroelastic loads. 

Analytically determined vibration modes of launch vehicle 
structures have been verified with full-scale vibration surveys 
performed in a manner similar to that developed for aircraft. With 
the increase in overall dimensions, this approach has become more and 
more expensive. Current major United States launch vehicles range 
from approximately 120 feet to 320 feet and the suspension of these 
vehicles to simulate the free-free end conditions of flight is, 
in itself, a major task. The simulation of configuration changes 
and various propellant conditions present additional severe restric- 
tions. The replacement of the full-scale vibration survey with an 
equivalent dynamic scale model test, can offer multi-million dollar 
savings in materials, manufacturing time and total test program 
time. 


One of the first major launch vehicle programs in the United 
States where the full-scale vibration survey was replaced by a 
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model test was the Titan III program funded by the U. S. Air Force. 
Several Saturn class dynamic models have since been examined, and it is 
probable that most, if not all, of the large launch vehicles of the 
future will be vibration tested by the use of dynamic scale models. 
Recent research performed at the Langley Research Center of the 
National Aeronautics and Space Administration indicates the definite 
feasibility and accuracy of this approach. 

Model testing is utilized exclusively for aeroelastic tests 
but model size is restricted by available wind tunnels. 

The purpose of this study is to develop modeling techniques 
that can be used for launch vehicle and space structure dynamic 
simulation when the parameters which define similarity between 
prototype and model cannot, due to physical or manufacturing 
limitations, be fully satisfied. In practice, all but the largest 
models have to make compromises in structural simulation and even 
very large models, using direct replica scaling, cannot duplicate 
all essential parameters. 

During the study it became obvious that there is need for a 
comprehensive summary of the scaling techniques applicable to the 
design of launch vehicle and space structure models. This report, 
provides that summary; describes general model scaling methods; 
developes basic similarity laws for various structures including 
liquid propellants and shell structures. A survey of material and 
manufacturing limitations and a summary of the economic aspects 
of model construction and testing is included. 
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2* REVIEW OF THE STATE OF THE ART 


From the number of problems which have been solved using 
scale modeling techniques, the potential and application of this 
useful method are obvious* Static and dynamic models are being 
used to study structural response of space vehicles ranging from 
transportation occurring during ground handling to docking maneuvers 
in orbit to lunar and planetary landing* 

A colloquim on the use of models and scaling in shock and 
vibration sponsored by the ASME in 1963 treated recent develop- 
ments. These included wave simulation, aerospace vehicles, 
materials and elastic and plastic response to transient loading* 

A symposium on aeroelastic and dynamic modeling technology was 
sponsored in 1963 by the Air Force Flight Dynamics Laboratory of 
the Research and Technology Division, Air Force Systems Command in 
association with the Dynamics and Aeroelasticity panel of the 
Aerospace Industries Association. This meeting included papers 
on such topics as transonic buffeting, ground winds, dissimilar 
materials, thermal scaling, support systems and excitation tech- 
niques, fuel slosh, fatigue, photoelastic techniques, plus 
applications to space vehicles and aircraft* 

All of the work on dynamic response of space vehicles is 
aimed at answering the same questions What are the dynamic 
stresses which we must design the vehicle to withstand? The 
current scale modeling technology for studying the various types 
of dynamic response is oriented toward those aspects for which 
analytical models are inadequate. The degree of reliance on the 
results of scale-model tests often depends on the willingness of 
the responsible engineer to rely on the results of analysis* 

This "willingness 11 is a function of schedule, available funds, 
past experience on the problem at hand, program penalties incurred 
if he believes analysis and program penalties incurred if analysis 
is wrong. 

The philosophy of the scale-model test to be run also depends 
on this "willingness". By natural cause and effect, a large 
amount of money can be spent to solve a very serious problem 
(e*g., the wind-induced oscillations of the Saturn V vehicle or 
the POGO oscillations of the Titan II-Gemini) . In this context 
"solve* 1 means "make sure that the mission does not fail". The 
associated requirement on the scale-model test is that, it must 
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fill the gap between the lack of confidence in analysis and the 
"solution of the problem". Obviously, the range of requirements 
is extremely broad. All scale model tests have one thing in 
common, however, everyone of them is run to provide the "missing 
information" required to "solve the problem." Since the avail- 
able knowledge of the problem area and the confidence in the 
results of analysis are usually closely correlated, the degree 
of similitude which is required in the scale-model test increases 
as the analytical model deteriorates. 

Three general types of models have been used for predicting 
dynamic response of space vehicles. One type is the obvious 
geometrically scaled model made of the same material as the 
prototype. The second is again geometrically scaled but of a 
different material; and the third is a distorted model whereby 
one or more of the similitude requirements is not preserved. The 
first type of model has been used successfully for tests to 
determine static and dynamic stresses, displacements, and strains. 
However, the extremely light weight structure of aerospace vehicles 
sometimes makes geometric scaling impossible. For example, the 
tank wall thickness on one launch vehicle is only 0.005 inch 
full scale. On a 10% scale model geometric scaling would require 
a wall thickness of 0.0005 inch. 

Scale models to predict dynamic response, therefore, often 
require departure from geometric scaling and the use of the 
different materials for model and prototype. 

The structural response problem (static, dynamic, rigid body, 
vibrational, buckling, linear, nonlinear, elastic, plastic) can 
be formulated in terms of the geometry of the structure (dimensions) , 
material properties of the structure (linear modulus, nonlinear 
modulus, damping, strain rate effects, density, Poisson's ratio), 
and forcing function applied to the structure (point loads, and/ 
or pressure as functions of time). Four basic physical dimen- 
sions describe all of the above. Dimensional analysis formulates 
the problem of gaining full scale information from models in terms 
of the similitude requirements which when satisfied will result 
in reliable scaling laws. 

It should be noted that the validity of a scale model test 
is measured by whether it provides the information required to 
solve the problem. The improvements in the state-of-the-art for 
scale-model tests can thus be accomplished in either of two ways: 

(1) Development of more sophisticated methods of building 
scale models and improving the simulation of the 
environments ; or 
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(2) Improvement of the methods of analysis and understanding 
of the physics of the problem, particularly in those 
areas which cause the most trouble in the design of 
scale-model tests. 

Dynamic models of launch vehicles have had widespread application 
to two major areas; vibration surveys of the natural modes of 
oscillation and aero elastic loads and stability testing. 


The most common method of dynamic model design and construc- 
tion is direct geometric scaling using similar materials for 
model and prototype (replica scaling). All prototype dimensions 
are scaled directly through the model-to-prototype characteristic 
length scale factor, n = x^/A^ » and prototype materials are used 

in model construction. 


In the design of launch vehicle dynamic models it is often 
necessary, due primarily to material or manufacturing limitations 
(Chapter 9), to employ limited distorted geometric scaling and/or 
dissimilar materials, limited distorted geometric scaling deviates 
from replica scaling in that the characteristic length scale factor is 
not maintained throughout the entire structure. Small scale models 
may require extensive use of these techniques to provide the required 
stiffness distributions and to ensure that the several model com- 
ponents are strong enough to withstand handling and assembly loads. 

In addition, distorted geometric scaling is almost a necessity for 
the design of aeroelastic models where the scale factor is dictated 
by the wind tunnel imposed limitations on model size. The require- 
ment of matching aerodynamic and structural scaling parameters at 
these small scale factors may dictate the use of dissimilar 
materials . 

Total geometric distortion is dictated when very large structures 
are to be simulated. In these cases the model may not even resemble 
the prototype but the overall dynamic characteristics are retained. 

2.1 Vibration Survey Models 

Some of the most significant vibration survey models have been 
associated with research programs performed at the Langley Research 
Center (LRC) of the National Aeronautics and Space Administration. 

Two models of the Apollo/Saturn V vehicle were constructed; one at 
1/40-scale (Fig. 2.1) and the other at 1/10-scale (Fig. 2.2). The 
1/40-scale model did not use direct geometric scaling and, in some 
areas, dissimilar materials were used. The main load carrying 
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Figure 2.2 Saturn V 1/10-Scale Model 
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structure of the 1/10-scale model booster stages was represented by 
direct geometric scaling using similar materials but direct scaling 
was not employed for the payload and launch escape system. Good 
agreement (Ref 2.1) between experimental and analytical lateral 
bending frequencies (Fig. 2.3a) and vibration modes (Fig. 2.3b) was 
obtained. A summary discussion of this model program is presented 
in Reference 2.1. 


— Simplified Pretest Analysis 

o 1/10-Scale Model 
□ 1/40-Scale Model 
— — Posttest Analysis 
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Figure 2.3 Lateral Vibration Characteristics, Saturn V Dynamic Models 


In addition to these models a 1/ 5-scale model of the Saturn I 
launch vehicle (Fig. 2.4) was constructed using direct replica 
scaling. 
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Figure 2.4 Saturn I Vibration Test Vehicles 
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Vibration test results for this model are compared with the 
results of a full-scale Saturn I vibration survey in Reference 2.2. 
First bending and cluster mode shapes and frequencies for the 
model and full-scale vehicle are presented in Figure 2,5. 
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Figure 2.5 Fundamental Modes at Maximum Dynamic Pressure-Weight, Saturn I 


Other launch vehicles that have undergone extensive analysis 
through the use of a dynamic scale model are the Titan IIIA and 
Titan IIIC. A complete ground vibration survey was performed at 
the Langley Research Center. The Titan IIIC model, in the test 
stand, is shown in Figure 2.6. 
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The model utilized direct geometric scaling with similar 
materials (replica scaling) and full-scale manufacturing tech- 
niques except that some deviations were allowed in the Transtage 
structure. The model was used to replace the usual full-scale 
vibration survey. A discussion of this test program is found in 
Reference 2,3. Comparison of representative model vibration 
modes is shown in Figure 2.7. Comparison of representative model 
frequencies is presented in Table 2,1. 

Table 2.1 Frequency Comparison, Titan IIIA 
and Titan IIIC 1/5-Scale Models 


Mode 

Titan IIIA 
5000 lb 

, Lift-Off 
. Payload 

Titan IIIC, Lift-Off 
26,000 lb. Payload 


Anal • f Hz • 

Expt. , Hz. 

Anal,, Hz. 

Expt., Hz. 

1st pitch 

11.7 

11.7 

5.9 

6.0 

2nd pitch 

23.9 

23.6 

13.7 

12.7 

3rd pitch 

30.7 

32.0 

18.2 

20.3 


Results of a statistical study conducted to establish the 
tolerances associated with the dynamic characteristics of the Titan 
IIIC are presented in References 2.13 and 2,14. Modal frequency 
tolerances for a 99% confidence level, are presented in Table 2.2. 


Table 2.2 Modal Frequency Tolerance, Titan III C 


Plane 

Time of Flight, Sec 

Values of Tolerances in % 

Mode 1 

Mode 2 

Mode 3 

Pitch 

0 

3.6 

7.3 

14.7 

Pitch 

105 

3.1 

5.3 

10.6 

Yaw 

0 

4.0 

7.5 

19.0 

Yaw 

105 

3.4 

5.6 

13.0 


Note that at modes as low as the third, a tolerance varying from 
10% to 19% is determined. 
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Figure 2.7 Mode Shapes, Titan IIIA and 1IIC l/5~Scale Models 
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Modal displacement and slope tolerances for the same configuration 
are presented in Table 2.3. 


Table 2.3 Modal Displacement and Slope Tolerance, Titan III C 


Description of Group 

Values 
Mode 1 

of Tolerances in % 
Mode 2 Mode 3 

All payloads, all flight times 

$ 

6 

21 

28 

pitch and yaw planes 

e 

36 

50 

45 

All payloads, all flight times 

<}> 

7 

12 

13 

pitch plane 

0 

47 

20 

12 

All payloads, all flight times 

4> 

6 

29 

121 

yaw plane 

0 

25 

99 

141 

5,000 lb payload, all flight times 

<f> 

3 

10 

9 

pitch and yaw planes 

0 

21 

37 

36 

26,000 lb payload, all flight times 

<P 

7 

51 

356 

pitch and yaw planes 

0 

48 

31 

39 

45,000 lb payload, all flight times 

$ 

11 

50 

227 

pitch and yaw planes 

0 

54 

147 

37 

All payloads, flight time = 0 sec. 

<i> 

8 

34 

19 

pitch and yaw planes 

0 

35 

47 

80 

All payloads, flight time = 53 sec. 

<P 

7 

18 

20 

pitch and yaw planes 

0 

25 

46 

17 

All payloads, flight time = 105 sec. 

4> 

2 

14 

36 

pitch and yaw planes 

0 

38 

35 

10 

Note 1, <j> = modal displacement, 0 = modal slope 

2 . Payload weights are full scale values . 


These tolerances represent average values of tolerances evaluated 
at individual sensor locations. 
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An analytical technique that can be used to improve the quality 
of experimental vibration modes was developed in conjunction with 
the Titan III Transtage test. Several modes with very little 
frequency separation resulted due to the complexity of the structure. 
These experimental modes were used as "assumed" modes for a vibration 
analysis and the coupled Transtage mass matrix was determined under 
the assumption that only inertial coupling exists. The uncoupled 
modes, computed by the vibration analysis, were then used as the 
Transtage modes for all further analyses. The uncoupled modes 
differed only slightly from the test modes but these differences 
were sufficient to provide modal orthogonality with respect to the 
mass matrix. 

Modal damping data (Ref. 2.14) obtained during the Titan IIIA 
and Titan IIIC 1/5-scale model test are presented in Figure 2.8. 



Figure 2.8 Modal Damping Data, -Titan IIIA 
and Titan IIIC 1/5-Scale Models 
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2.2 Aero elastic Models 

Aeroelastic models of launch vehicles have been primarily 
constructed to investigate the response, and imposed loads, of 
vehicles subjected to ground winds and to random transonic buffet 
excitation. 


2.2.1 Buffet Models 

Available wind tunnels impose a size restriction on transonic 
buffet models and typical models are restricted to lengths of 10 feet 
to 16 feet. Figure 2.10 shows a 7% scale aeroelastic buffet model 
of the Titan III vehicle with various payloads . This model was 
scaled in accordance with the scale factors dictated by direct 
geometric scaling using similar materials but was not constructed 
with replica scaling. Only the simulation of bending modes was 
necessary to determine buffet response. The bending and shear 
stiffnesses were simulated with a machined metal tube of smaller 
diameter and higher wall thickness than a directly scaled replica. 

The model is shown in Figure 2.9. A basket-like machined pattern 
of cut-outs provided the proper ratio of bending to shear stiffness. 



Figure 2.9 Titan III 77» Buffet Model Construction 
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C. 260-in. Bulbous Payload 


Figure 2.10 Titan III Aeroelastic Buffet Model 
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The more important scaling parameters (Ref 2.4) of this Titan III 
model are summarized in Table 2.4. The subscripts m and p refer to 
model and prototype, respectively. 

Table 2.4 Scaling Parameters, Titan III 7% Buffet Model 


Variable 

Scaling Relation 

Length 

£ /£ » 0.07 = n 

m p 

Total Mass 

M /M = (0.07) J 
m p 

Velocity 

V /V = 1 
m p 

Frequency 

w /w = 1/0.07 
tn p 


The test operating medium and conditions were such that the 
relations between the non-dimensional parameters Strouhal Number 
(oj£/V) and Reynolds Number ( pV £/ U) were 


S/S = (w/w )(£/£) (V/V) 
m p mpmppm 


Re /Re = (p / p ) (V /V ) ( £ /. £ ) (U /u ) 
mp mpmpmppin 


n 


A complete discussion of the model and test can be found in 
Reference 2.4. Comparison between model and scaled full-scale 
bending frequencies is shown in Table 2.5. 


Table 2.5 Bending Frequencies, Titan IIIC 



Pitch, Hz. 

Yaw, Hz. 

Model 

Full-Scale 

(scaled) 

Model 

Full-Scale 

(scaled) 

1st Bending 

24.3 

23.2 

23.3 

20.5 

2nd Bending 

40.0 

40 .3 

49.1 

39.7 

3rd Bending 

76.8 

71.8 

150.2 

79.7 


Typical test results are presented in Figure 2.11. 
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Non "Dimen s iona 1 
Moment At Stn. 

Normalized to 
Maximum Moment 
in Titan IIIA 
Configuration 



Mach Number 


Figure 2.11 Normalized rms Oscillatory Core Bending Moment, 
Titan III Model 


Aeroelastic models of the Saturn I launch vehicle have been 
built at 2% and 8% scale using dynamically similar construction. 
Complete structural replica scaling was not attempted. Test 
results for these models are presented in References 2.5 and 2.6. 
The 2% model was used primarily for the measurement of aerodynamic 
damping but it can be regarded as being characteristic of present- 
day buffet models. 

The scaling approach used to transform model test results into 
usable full-scale values is developed in Reference 2.7. This 
document also presents test results for a simple buffet model (with 
two payload shapes) characteristic of single-body launch vehicles. 
References 2.8 through 2.10 present recent test results obtained 
during an extension to the Titan III buffet test program. Several 
configurations and payloads are considered. 
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2.2.2 Wind-Induced Oscillation Models 


Several launch vehicle models have been built to predict dynamic 
loads induced by ground winds . Ground wind induced oscillations are 
associated with periodic vortex shedding on the vehicle and umbilical 
tower. Simulation, or at least resonable approximation, of the full- 
scale Reynolds Number is an important requirement and usually 
necessitates testing in media other than air (e.g. , Freon). Proper 
simulation of the test gas-model density ratio then requires higher 
values of model mass than are dictated by direct density scaling. 

One of the most significant ground wind models was the 7.5% 
aeroelastic Titan IIlA model and launch tower shown in Figure 2.12. 
The model was mounted on a turntable to permit variation of the 
relative wind direction. 

As it was desirable to match full-scale Reynolds Numbers as 
closely as possible, the model was tested (using Freon 12 as the test 
medium) in the Langley Research Center 16 foot variable pressure 
tunnel. This use of Freon as the test medium imposed certain 
restrictions. The imposed density and viscosity ratios were 


P / P = 3.7 and u /u = 0.757. 
m p m p 

Even with Freon it was impossible to match full scale Reynolds 
Numbers using practical materials for model construction without 
violating mass or stiffness scaling parameters. The test Reynolds 
Numbers were 0.45 of the full-scale values. Some other important 
scaling parameters for this Titan III model are presented in 
Table 2.6 


Table 2.6 Scaling Parameters, Titan III 7.5% 

Aeroelastic Wind- Induced Oscillations Model 


Var iable 

Scaling Relation 

Length 

£ /£ = 0.075 = n 

m p 

Velocity 

V /V = 1.228 
m p 

Frequency 

co /to =16.37 

m p 

Dynamic Pressure 

Vs - 5 - 576 

Density 

Vs = 3 - 7 
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Figure 2.12 Titan IIIA 7.5% Aeroelastic Model with Umbilical Mast 
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Note that the Strouhal Number (w£/v) was scaled exactly. The mass 
ratio was 


M /M 
m p 


(.u IJ 3 (PJ pj = 1*55 X 10" 3 . 
m p m p 


These scaling parameters, with the exception of the Reynolds 
and Mach Numbers, represent a correct set for aeroelastic scaling. 

The test Mach Numbers did not equal the full-scale Mach Number but 
this discrepancy was shown to be negligible through calibration 
tests in which the Reynolds and Strouhal Numbers were held constant 
while the Mach Number was varied. 

A summary discussion of this test program is presented in 
Reference 2.4 and a complete report is to be found in Reference 2.11. 
Response of the vehicle or umbilical structure (if affected by vortex 
shedding) is primarily a function of the Strouhal Number. Figure 2.13, 
taken from this reference, shows the dimensionless acceleration 
power spectral density measured for typical Titan III payloads and 
umbilical mast configurations. 



Figure 2.13 Normalized Distribution, Titan III Wind-Induced 
Oscillations Model 
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Figure 2.14 gives an indication of the effect of 
direction on the dynamic bending moment response 
umbilical mast at the Strouhal Number of maximum 


relative wind 
of the Titan III 
dynamic response. 


Relative Wind 



Relative Wind Azimuth ((3) 

Figure 2.14 Contours of Constant Dynamic Moment Coefficient, 
Titan III 7.5% WIO Model 


The dynamic bending moment coefficients used in this chart were 
determined from the relationship 


V - M dyn /qLR 
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where 


= measure ^ dynamic bending 
q - tunnel dynamic pressure. 


moment. 



FL 


L = length of mast, L 


and 


R 


J' <|>(x)W(x)dx 


where x is the mast coordinate. The function <j)(x) represents 
the normalized first bending mode shape (which accounted for 
practically all of the response) and W(x) is the width of the 
mast perpendicular to the wind direction. The parameter D 
(Figure 2.14) used in computing the reduced velocity (inverse 
Strouhal Number) is the width of the tip of the mast. 

Other aeroelastic models designed to investigate the effects 
of ground winds (Fig. 2.15) have included a 3%-scale Saturn V 
model, a 15%-scale model of the Scout vehicle and a 20%-scale 
Jupiter vehicle. The use of the Langley Research Center 16-foot 
variable density transonic tunnel permitted close simulation of 
the full-scale Reynolds Number except for the Saturn model. In 
this case the model was tested at approximately one-third of the 
full-scale Reynolds Number (Ref 2.2). 

The maximum lateral oscillatory bending moment as a function 
of the inverse Strouhal Number (Ref 2.2) for the 3%-scale Saturn V 
model is shown in Figure 2.16. A large response over a limited 
range of inverse Strouhal Number is indicated for the lowest 
damping case. This peak is seen to vanish as damping is increased. 

A summary of publications dealing with the effects of ground 
winds on launch vehicles will be found in Reference 2.12. 
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Figure 2,16 Effect of Damping, 3% Saturn V Model 
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3. DYNAMIC MODEL SCALING PHILOSOPHY 


The basic requirement governing the design of any dynamic 
model is that of similitude between the model and the full-scale 
structure (or phenomenon) which is to be simulated. This require- 
ment must be satisif ied for any of the time-dependent physical 
processes or characteristics that are of interest to the experi- 
menter. The analytical tool generally used to develop the 
similitude laws governing model design is dimensional analysis . 

The development of this method (Ref 3.1, 3.2, 3.3, and 3.4) is 
attributed to Buckingham and Rayleigh. Dimensional analysis is 
based upon the hypothesis that the solution of a physical problem 
is expressible in the form of a dimensionally homogeneous equation 
written in terms of the variables affecting the problem. Once 
these variables are specified, a group of dimensionless scaling 
parameters can be derived. Similitude is achieved when the scaling 
parameters are identical for both the model and the prototype. 

3 . 1 Scaling Laws for Structural Vibrations of Launch Vehicles 

The prediction of the vibrational characteristics of launch 
vehicles through the use of scale models is based upon the assump- 
tion that both model and prototype are governed by the same physical 
laws. From these physical laws can be established the modeling laws 
which mathematically relate the corresponding characteristics of the 
two systems. The non-dimensional form of the system governing 
equations must be the same for both model and prototype. 

In general, structural models can be separated into three 
classes; the particular model choice is governed by the complexity 
of the prototype structure, the type of exciting forces and the 
frequency range to be investigated. A graphical representation of 
the classes is shown in Figure 3.1. 

If the exciting frequencies are below the first few structural 
frequencies many of the details present in the prototype vehicle 
can be neglected, the structure can be represented by mass-spring- 
damper elements and, within this limited frequency range, the model 
will exhibit faithful dynamic similarity. As the exciting frequency 
spectrum increases a more complete geometrical similarity is 
required as the prototype vehicle no longer behaves as a single 
unit but rather becomes a multi-branched structure. 
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Figure 3.1 Dynamic Model Classes 


The form of applied excitation also governs the choice of 
model. If the applied forces are approximately uniform over the 
surface, geometric scaling of skin panels is a reasonable approach. 
For propagating excitation (e.g. , noise) it becomes necessary to 
accurately scale geometry, mass and stiffness characteristics. 

There exist two methods of predicting the behavior of physical 
systems. One is by means of analysis or computer simulation whereby 
a mathematical model or an analog is created of the physical system 
in question. The second method is by tests on a prototype or a 
properly scaled model. Considering the model test, the problem is 
to identify the physical parameters which govern the phenomenon of 
interest, to cast them into groups of non-dimensional variables, 
to perform tests on properly scaled models, and to apply the results 
from these tests to prediction of the behavior of the prototype. 
Limitations exist for both systems, Assumptions made in the 
mathematical model may be too restrictive or inaccurate, thereby 
limiting the range of applicability. More complete mathematical 
models may be too complex for economical solutions. Scaled models 
may not be satisfactory if all the important variables are not 
identified and/or if proper similitude is not provided. The theory 
of dimensional analysis is based on the Buckingham Pi theorem from 
which scaling laws are evolved. 
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The measure of physical phenomena can be expressed in units 
of the independent quantities of mass, length, time, and tempera- 
ture either singly or in multiple combinations.* Assume a physical 
phenomenon y depends on several independent variables, x., • • • 

x in some unknown manner. If all the independent variables are 
known, it is possible to predict the outcome y for a given set 
x- , x^» . . • x n by means of a properly scaled model test, even 
if the functional relationship between y and its independent 
variables in unknown. 

The following theory underlies the derivation of the scaling 
laws. The functional relationship between y and all the variables 
on which it depends can be expressed in general as 

y = f(x x , x 2 , . . . x n > (1) 


If the exact nature of the function f were known, then this 
would be the mathematical expression of the physical law expressing 
the dependence of y on x^, x^, . • * x^. This physical law is 
independent of the units used in the measurement, and applies 
just as well to the model as it does to the prototype. Using 
the subscripts m to apply to the model and p to apply to the 
prototype, we have the relationships 

y » " £(x l»' x 2m* ' 

y p ' £ °V x 2p* • 


. . X ) 

nm 


. .. x ) 
np 


( 2 ) 


Since the variables in a physical problem can be expressed 
in terms of the four basic dimensions of mass, length, time, 
and temperature, then the n + 1 variables (y , x^, x^, . . • x ) 
can be combined into (n + l)-4 = n-3 dimensionless groups. a 

Remembering that the physical law is the same for model and 
prototype. Equation 2 can then be rewritten 


Tr- 

im 

Tr- 

ip 


F(rr 2m , • * * Tr(n-3)m) 

F(Tr 2p» 7r 3p’ * * * Trl ^ n “ 3 )p) 


( 3 ) 


where the dimensionless groups TT, , 7t , contain the dependent 
variable y, and the other dimens ionlesl groups contain the 
independent variables. If the model is so constructed that the 
similitude re qui r ements 

*Electrical and magnetic phenomena require additional units. 
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7r(n-3)m = ir(n-3)p 


are satisfied, then it must follow from Equation 3 that the scaling 
law for the dependent variable is 


7T. 


lm 


TT- 


IP 


(5) 


Since y (contained in TT ) 
test, the dimensionless quantity 
determined experimentally. 


can be measured from the model 

tt * tl can therefore be 
Ip lm 


Thus knowing tt^ , the numerical value of the physical quantity 
y contained in tt Ban be calculated. Equation 5 thus gives the 
dlsired scale factBr for the dependent variable y^. 

With four basic dimensions in the problem (mass, length, time, 
and temperature) , four scale factors may be chosen arbitrarily to 
suit test requirements. The remaining scale factors will then be 
expressed in terms of one or more of these four. 


In general, if all the variables on which a particular 
phenomenon depends are known, dimensional analysis will give the 
desired scaling law for the dependent variable and all the 
similitude requirements necessary to ensure its theoretical 
validity. However, two important practical difficulties usually 
must be overcome to develop a workable scaling law* The first 
difficulty is the practical impossibility of satisfying all the 
similitude requirements that are theoretically necessary (Eq. 4) . 

An approximately correct scaling law, however, can still be obtained* 
Only by experimentation can the importance of the similitude 
requirements that are necessarily neglected, be determined. In 
any case, empirical corrections to the scaling law can usually 
be developed. For example, if a particular similitude require- 
ment, tt « tt , proves to be important, and yet cannot be 
satisfied, the p dif f iculty can be overcome by determining experi- 
mentally how tt^ varies with tt^, keeping all other TT terms constant. 

The second difficulty is that a complete knowledge of all 
the variables may not exist. Thus, the unsatisfied similitude 
requirement for each unknown variable may cause deviations from 
the scaling law. 
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The study of the modeling parameters is applicable to the 
determination of lateral vibrations of aerospace vehicles* vehicle 
stages* and orbiting space structures such as antennas* mirrors 
and space stations * The application of scaling techniques has 
been used to study a wide variety of mechanics problems. The 
work in this study investigates* in a quantitative manner* some 
applications of unconventional modeling techniques. The purpose 
of this approach is to identify and to evaluate those techniques 
that can be used on structures and structural problems that do 
not lend themselves to the more common modeling techniques. 

The formulation of a problem in dynamics response requires 
synthesis of three related parts. These are 

1. Geometry of the structure and its components. 

2. Mechanical properties of the material and characteristics 
of joining techniques used in the structure. 

3. Loads and environment of the structure. 

Dimensional analysis is used to identify the numerical values 
of scaled indepdendent variables (geometry and materials properties) 
required to provide complete similitude for lateral vibrations of 
space vehicle. Requirements of model geometry and model material 
cannot be satisfied using conventional modeling techniques. It 
is the intent of this study to determine the values of materials 
properties (such as density and modulus) required for the preserva- 
tion of lateral vibration similitude. The limits imposed on model 
size are identified by the availability and economy of unusual 
techniques . 

A study of the convenience of the four arbitrary scale 
factors can be undertaken to identify which combinations are most 
convenient for the particular aspect of the problem being modeled. 
One requirement for the set of properties (independent variables) 
chosen is that the basic dimensions of M* L, T* and 0 must be 
obtainable by combinations of those properties, A very common 
choice is as follows! 

Symbol Property Dimension 

H Length L 


p Mass Density ML 


Scale Factor 


m 


m 


= n„ 


m 


M L~ 3 
m m 

M L~ 3 p 
P 
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Symbol Property 

v Velocity 


Dimens ion 



e 


Temperature o^ 


Scale Factor 



L T 
m 


-1 


m 


L T 
P 


-1 


= n„ 



n 4 


Once chosen, these factors must be used for scaling on all of the 
aspects of the given problems 

Instead of the list given above, any set of independent 
variables which contains all four of the basic dimensions can be 
chosen. For example, one can choose a pressure variable instead 
of mass density since Pv^ has units of pressure. In the end, 
however, the similitude requirements for all the other independent 
variables and the scaling laws for the dependent variables are 
expressed in terms of the four scale factors n^, n^, n~, and n^ 
which have been chosen. Thus, this choice determines the 
specific difficulties which will be encountered in building the 
scale-model, A particular difficulty can be exchanged for a 
less-formidable difficulty by altering the choice of arbitrary 
scale factors. 


Candidate fabrication techniques for metal models ares 
electroplating copper or mickel on a wax mandrel (a system used 
for very thin shell specimens for buckling experiments); vacuum 
deposition of thin films on lost wax and chem-milling (currently 
used on aluminum and titanium aerospace structures). 

Application of nonmetallic materials (plastics, fiberglass 
composites, and sandwich composites) for models are possible. 

It is important to identify those materials which can provide 
appreciable variation in other properties such as modulus and 
density and still maintain the required geometric similitude. 


Producing models of materials different from the prototype 
generally means a difference in material paramters (E, p, and v) . 
For example, let us consider the following similitude requirement. 


2 2 

P l p it 

m m P P 

2 2 ’ 

E t Z E t' 

mm P P 


D E .N V t R DIVISION 



MCR-68-87 


3-7 


Assume that all times for model and prototype are required to be 
the same for a 10% model. We then have 


1 = 



1 

100 


Obviously, this cannot be satisfied using the same material for 
the prototype and model. If a broad choice of material is avail- 
able, we can choose for the model 


£ 

m 


E 

p , a very flexible material 
10 


m 


10 p 

P 


a very dense material 


A phenomenon that can be important in modeling using plastics 
is the presence of viscoelastic behavior. If the model material 
exhibits a strain-rate effect during vibration studies, but the 
prototype material does not, then the modeling laws must account 
for this violation of the similitude requirements. 


Requirements for adequate scaling of joint details (including 
friction, roughness, prestress, and internal forces of fasteners) 
is a formidable task. The need to make small-scale exact replicas 
of structural joints results from the lack of an analytical or 
empirical model of the dynamic behavior of these joints. 


3 . 2 Basic Dimensions for Scaling Analysis 

For the several classes of problems to be considered, we 
will first establish three arbitrary scale factors from which 
we can isolate mass (or force), length, and time. Although the 
primary interest is dynamic scaling, we shall not neglect the 
static problem. The procedure is to build the dynamic problem 
from the static and to have forced vibrations related to free 
vibrations. 

To indicate many of the terms that should be considered in 
structural scaling, Table 3-1 is presented in dimensional form. 
Both the engineering dimensions of force (F) , length (L) , and 
time (T) are given as well as the basic dimensions of mass (M) , 
length (L) , and time (T) . Converting from one system to the 
other we use for the definition of force mass times acceleration 
or symbolically 

F = MLT -2 


TVO/V 
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Table 3,1 Physical Quantities and Dimensions 


Quantity 

Symbol 

Engineering Units 

Basic Units 

Length (Any linear dimension) 

a 

in. 

L 

L 

Radius 

R 

in. 

L 

L 

Thickness 

h 

in. 

L 

L 

Displacement 

w 

in. 

L 

L 

Dynamic Response Amplitude 

X 

in. 

L 

L 

Area 

A 

in *? 

L 2 

L 2 

Volume 

V 

in. 3 

L 3 

L | 

Moment of Inertia 

I 

in.^ 

L 4 


Polar Moment 

J 

in *4 

L 4 

L 4 2 

Force 

p 

lb. 

F 

MLT 

Moment 

M 

in. lb . 

FL 

ML 2 T“ 2 

Pressure 

q 

lb / in 3 

FL-2 

ml - 1 t -2 

Stress 

a 

lb /in 2 

FL “2 

ml-i t - 2 

Modulus of Elasticity 

E 

lb /in 2 

FL” 2 

ML” 1 !” 2 

Shear Modulus 

G 

lb /in 2 

FL -2 

ml - 1 t -2 

Unit Mass Density 

P 

lb ,sec 2 /in 4 

FL“ 4 T 2 

ML“ 3 

Unit Weight Density 

Y 

lb /in 3 

FL-3 

ML” 2 !” 2 

Velocity 

V 

in/ sec 

LT-1 

LT" 1 

Acceleration 

a 

in/sec2 

LT~ Z 

LT -2 

Gravity constant 

g 

in/ sec 2 

LT“ 2 

LT”” 2 

Time 

t 

sec 

T 


Frequency 

CO 

1/sec 

T“1 

T“1 

Beam Extension Rigidity 

AE 

lb 

F 

MLT" 2 

Beam Bending Rigidity 

El 

in2lb 

FL 2 

ML 3 T“ 2 

Beam Torsion Rigidity 

GJ 

in2lb 

FL 2 

ml 3 t“ 2 

Beam Shear Rigidity 

kAG 

lb 

F 

mlt“ z 

Shell Extension Rigidity 

K 

lb /in 

FL -1 

MT 

Shell Flexural Rigidity 

D 

in. lb. 

FL 

ML 2 T“2 

Poisson Ratio 

V 

— 

— 


Wave Number 

n 

— 



— 

Coordinate Angle 

$ 

— 

— 

— 

Arbitrary Scale Factor, n^ i 
Dimensionless Quantity, 7T^ 

Subscript m (model) 

Subscript p (prototype) 

- 1, 2, 
i - l, : 

... 

- f i . . 
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Several classes of problems become apparent when considering 
launch vehicle structures. They may be separated into beam-type 
and shell-type structures, free vibrations or forced vibrations, with 
liquid or solid propellants, various boundary conditions and several 
gravity-type loadings. Some of the previous general discussions 
regarding tt terms and arbitrary scale factors are considered in 
more detail for several kinds of structures, loadings and responses. 

To isolate the basic dimensions using arbitrary scale factors 
n^, n^ and n_ we will establish, the scaling ratios of F, L, and T 
for model and prototype. The usual arbitrary scale factor is the 
length factor and will be designated n^ so that 

i l 


m m 



This can be written for the model quantity containing as 

L = n.L (6) 

m 1 p ' 


Using the ratio of modulus of elasticity as another arbitrary 
scale factor we have 


n„ 


m 


F L 
m m 


-2 


E r t " 2 

p F L 
P P 


From Equation 6 for n^ this can be expressed as 


n 2 F n l 


jn -2 
F 

P 

and solving for the model force 


F * n. n 0 F 
m 1 2 p 


( 7 ) 


To incorporate the time scale factor the third arbitrary 
scale factor will be ratio of mass densities so that 


m 


F L “ 4 T 2 
mm m 

F L ” 4 T 2 
P P P 
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By using Equations 6 and 7 we can express n^ as 


2 -4 T 

n 3 - (n x n 2 ) ) jn_ 

T 2 
P 


and solving for T we have 


r . „ „ “ 1/ 2 1/2 

X « n. n„ n_ T 

m 12 3 p 


( 8 ) 


Based on the previous discussion of ir terms it is now possible 
to scale several kinds of problems, and to design properties into 
the model depending on the quantity to be measured. It is important 
to note that independent variables are scaled, dependent variables 
measured, and this measured value can be correlated to the prototype 
by the scaling laws. Important pieces of information available 
from the it terms include the ranges of physical values that must 
be provided or may be impossible to satisfy. 

3,2.1. Static Beam. Scaled Geometrically 

For the static problem either stress or deflection are the 
usual dependent variables of interest. We can isolate the 
important independent parameters and write the following: 


f(p, 

pg, M, SL, E, v) 

(9a) 

g(p* 

pg, M, l, E, v) 

(9b) 


Dependent ir terms 



Independent IT terms 



MARTIN MARIETTA COR RORATION 

DENVER DIVISION 



MCR-68-87 


3-11 


*5 



TT 


6 



1 


Applying the scaling laws ir^ m * 

Model dimensions from the it terms in terms of arbitrary factors 
become for the dependent variables 


u 

m 


m 



n.u 
1 P 


n 2 °p 


and for the independent variables 


P 

m 


Pr-8, 


mm 


M 

m 


E l 
m m 


E 

„ 2 

P 

SL n 


P 

E 

a 

IQ 

Z2. 

E 

1 

P 

m 

E 

i 3 

m 

m 

E 

P 

n 3 

p 


“lVp 


P«*r 


-1 

n, n 


2°P S P 


M » n, n_M 


V = V 
m p 


(10a) 

(10b) 

(Ha) 

(lib) 

(11c) 

did) 


Model deflections scale directly as the length scale factor and 
stresses as the modulus scale factor providing the other quantities 
are properly scaled. 

If the weight of the structure is not important then can 
be neglected. If, however, the structural weight is important 
similitude requires that tt^ be satisfied. 

3.2.2 Static Beam, Not Scaled Geometrically 

If for some purpose of convenience the model does not faith- 
fully reproduce the prototype in a geometric way then other 
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properties need to be specified and scaled. Geometric scaling 
automatically provides similitude for the additional terms 
involving properties of areas. 


u “ f(P, pg, M, Si, E, V, El, GJ, kAG) 

C " g(P, pg, M, Si, E, v, El, GJ, KAG) 
Dependent tt terms are 


tt. 


u 

7 


TT, 


£ 

E 


Independent tt terms become 

TT 

it 3 2 

EJ i 


TT, 


P8 £ 


TT, 


M 

EE 3 

P 


TT-, “ 


TT, 


El 

Pf, 2 

GJ 

?i 2 

kAG 

P 


(10a) 

(10b) 


For tt - tt and applying this to the above it terms we can scale 
model^quan^ities for the dependent variables from Equation 10a 
and 10b 
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u 

* n- 

u 

m 

1 

P 

cr 

88 n 0 

a 

m 

2 

P 


and for the dependent variables using Equations 11a through lid 
and including tt^, TTg and we have for the model 


P ■ n. n„ P 
m 1 2 p 


m m 


n l' 1 n 2 Pp *p 


H 

3 

n, n„ M 

m 

1 2 p 

v 

V 

m 

P 

p a 2 

E I 


m m 

P Z 2 

P P 

P l 2 

G J 

in m 

m m 

P 0 2 

P *P 


p p 


n* n_E I 

1 2 p p 


G J - n. n„G J 
P P 1 2 p p 


(12a) 


(12b) 


k A G 
m m m 


~ k A G 
Pp P P P 


n, n„k A G 

1 2 p p p 


(12c) 


Quantities from the geometrically scaled model are retained. 

In addition the bending, torsional, and shear rigidities are included 
to ensure their proper scaling. If, however, one or more of these 
additional terms is not considered important, depending on the 
loading condition, then that scaling can be disregarded. 

3.2.3 Free Vibrations. Geometrically Scaled Beam 

Another class of problems involving beams and beam-like 
structures is free vibrations. Many launch vehicle structures 
are large beams of circular cross section. Boundary conditions, 
and loads change on the beam from prior to launch, to liftoff, 
to staging, etc. 

The usual depedent variable for free vibrations are the 
natural frequencies uu If, however, stress or displacement are 
also required, the previous treatments given by Equations 10a 
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and 10b apply, providing the important similitude requirements 
have been met. For free vibrations of a beam subjected to a load 
(loads) P and with unit mass density p the important independent 
variables can be expressed as 

CO - f (P, pg, JL, E, V ) 

The dependent n term for the natural frequency co is given by 


tt. 





In a manner similar to the static problem the other tt terms 
become 


TT 


2 



TT 


3 


Pg 


JL 

E 


n, 



1 . 


Making tt ■ 
for the model as 


Tt for dynamic similitude we find the frequency 


co 


m 



co 

m 





(13) 


This result is to be expected after examining Equation 8 for model 
time T m , The frequency has the dimensions of T”1 (referring to 
Table 3-1) and obviously Equation 13 is the reciprocal of Equation 
8. When modeling with identical materials ■ 1, n_ ■ 1 then 
the well known expression for model and prototype frequency 
relationship (l/n^) is apparent from Equation 13. However, as 
shown for the independent variables given below, when the mass of 
the structure compared to the applied load (loads) P is large and 
must be scaled then similitude is violated when n ^ = 1 and model 
and prototype are of the same material. 
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For the independent variables from Equation 11 


P 

m 



P 

P 


P-8, 


mm 


V 1 n 2 P p g p 


V m 


V_ 


If the external load (loads) P are significantly greater than 
the weight of the beam then scaling the model loads by n^ x 
prototype load (and n£ * 1) then model and prototype can be of the 
same material and frequencies scale as given. It is theorectically 
possible, by examining Equation 13 to select any frequency ratio 
that appears convenient. The n^ term is the length factor and 
can be manipulated rather easily. Both n£ and n^ refer to the 
material properties of a modeling substance. Manipulations such 
as centrifuges, composite materials, additional masses are some 
of the devices used to modify these arbitrary scale factors. 


3,2.4 Free Vibrations. Beam Not Scaled Geometrically 

By using the three additional terms introduced for beams 
not scaled geometrically and by using the procedure outlined above 
this class of problem can also be scaled. 


3.2.5 Forced Vibrations. Geometrically Scaled Beam 

For forced vibrations the dependent quantity of interest is 
generally the amplitude of vibration. Both stress and displacement 
have been discussed and modeling for these quantities has been 
established, ^We will use the hat symbol to denote the forcing 
frequency as ou and the peak value (sinusoidal excitation) of the 
forcing function as P. Other symbols refer to quantities defined 
in Table 3-1. 


x • f (P , pg , H/ E, v, U)f P) 

If dynamic similitude is provided then the dependent variable x 
can be nondimensionalized as 



To ensure dynamic similitude the independent variables when made 
nondimensional the corresponding terms must be made numerically 
equal. These variables are given by 
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P 



IT- * (Jdi P 
D E 

A 


P 



Previously developed equations permit us to write directly for the 
dependent variable, the amplitude of vibration 


m 


n- x 
1 P 


For tiie several independent terms 

„ 2 
P ■ n. n 0 P 
m 1 2 p 

-1 

p g ** n. n 0 p g 
K m b m 1 2 K p°p 

v = v 

m p 

"" -1 1/2 - 1/2 A 

w ■ n, n„ n 0 w 

m 1 2 3 p 


P ■ n. n„ P 
m 1 2 p 


(14a) 

(14b) 


To scale the amplitude of the response by the length scale 

factor the model must be vibrated as given by w with a peak force 

of P . If the applied model load (loads) meetsthe requirement of 

P and the weight of the structure is not significant then similar 

materials for model and prototype provide well defined scaling. If 

however, there are no loads, then tt^ must be satisfied and the term 

p E becomes important, 
mm 

3,2,6 Forced Vibrations. Beam Not Scaled Geometrically 


As discussed in 3,2.4, no new techniques or terms are required. 
Those rigidities not automatically satisfied by geometric scaling 
(bending, torsion, and shear) are scaled with the n^ and n^ terms to 
appropriate powers as developed by Equations 12a to 12c. 
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3.3 Matrix Approach to Compatible Scaling Parameters 

A generalized matrix approach to dimensional analysis can be 
used to develop scaling laws for vibrations of launch vehicle 
structures. The technique is straightforward anc * reduces the 
problem to one of matrix algebra. 

In principle, a physical problem will involve a functional 
relationship among n variables necessary to describe the phenomonon. 
For example 


v „ ■ £ <V v 2 Vl>> ‘ 15) 

or equivalently, 

g(V r V 2 , . . . , V n ) = 0 (16) 

Dimensional analysis presupposes the form of Equation 16 to be 
that of a product of the pertinent variables, exponentiated in 
accordance with their unknown contribution to the product. 
Theoretically, there are n functional relations of the type shown 
in Equation 15 (or 16) • However, all n of these relations will not 
be independent as will be shown. Further, if Equation 16 is 
dimensionally homogeneous (that is, it does not depend on the 
fundamental units of measurement) it can be reduced to a relation 
among a complete set of dimensionless products, or scaling 
parameters, *rr ^ > 

F(t[^> 7^2* * • • * TT^) = ^ (17) 

In essence, this is a statement of Buckingham’s Theorem (Ref 3.5). 

When a force, length, time system of measurement is used, the size, 

N, of a complete set of dimensionless products is equal to the number, 
n, of necessary variables minus the number, d, of fundamental dimen- 
sions (Ref 3.6) . Thus, there are N independent relations of the 
form of Equation 16. Equation 17 does not imply that the functional 
relationship is known, but merely that one exists. Calculation of 
the dimensionless products , tt . , is accomplished through dimensional 
analysis. However, the functional relationship, F, among these 
products must lie within the realm of experimentation. Whatever 
the form of Equationl7, it must be satisfied by both the model and 
the prototype. 

A dimensionless product can be expressed in the form 



i=l 


N (18) 
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A dimensional equation (one expressed in terms of fundamental 

dimensions only) can be obtained from Equation 18 if each variable 

V_^, is expressed in terms of its fundamental dimggsions. For 

example, if the acceleration of gravity is the n variable, 

thep V * - LT~^ where the * indicates dimensionality, 
n 


Let 

V i - A \ ^ CW) 

k=l 

th 

where D is the k fundamental dimension. The exponents, q. ., 
can be determined if Equation 15 (or 16) is expressed in 
product form 


n 


1=1 1 


- 1 


( 20 ) 


where K is some constant. 


.th 


This is equivalent t^the j dimensionless product of Equation 
4, where is the j column of q_^ . 

Substitution of Equation 19 into a dimensional form of 
Equation 20 yields 


n 


rr 

i=l 


' d P, . 

n 

k=l 


, c i 


= 1 


( 21 ) 


In order for Equation 21 to be satisfied, all exponents must 
vanish. Therefore, 


[P ki ] tc ± l - 0. 


( 22 ) 


Equation 22 represents d linear, homogeneous equations in n unknowns. 
For a well formulated analysis the rank of [P..] will be d, so 
that (n-d) - N independent solutions exist. Tnis is consistent 
with the number of dimensionless products in a complete set (see 
Equation 17) • Equation 22 can be solved by several techniques 
of matrix algebra, one of which is outlined as follows: 

1) Remove N columns of [P] by post-multiplication with an 
(n x d) transformation matrix [Z^] > forming a (d x d) 
matrix [P^], 

[P L ] = [P] [Z x l (23) 
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where [Z^] must be chosen so as not to render [P^] 
singular. The remaining columns of [Pj then form a 
(d x N) matrix [P^l so that 

[P 2 1 - [Pi [ Z 2 1 (24) 

where [Z^] is an (n x N) transformation matrix. 

2) Let [Ql be an arbitrary (N x N) matrix whose columns 
are mutually orthogonal (such as the (N x N) identity 
matrix). Then [Ql represents an arbitrary choice of the 
N independent q f s of Equation 18, forming N linearly 
independent vectors. A matrix equation analogous to 
Equation 22 can be written f or the exponents of the N 
dimensionless parameters 

[P x ] [U] + [P 2 3 [q] = o 

where the vectors of [U] are the unknown q. .Vs. 

Therefore, since [ P^] is non-singular, 

[U] = -CP^- 1 [P 2 ] L Ql • 

3) The exponents of the pertinent physical variables, V., 
can be determined by assembling the [U] and [Q] 
matrices according to 

[ q ± j ] = [Z 1 ] [U] + [Z 2 ] [Q]. (27) 

4) With the solution of Equation 27, a complete set of N 
dimensionless scaling parameters is determined. 

3.3.1 Example of Scaling Parameter Derivation 

Consider the vibration of a slender beam with bending and 
shear stiffness. The functional relationship between the physical 
variables is 

a) = f (m, £, AG, El) 

or f(m, £, AG, El, oj) = 0 

where w = natural frequency, 1/T 


(25) 


(26) 
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2 2 

m = mass/unit length, FT /L 
l = length, L 
AG = shear stiffness, F 

2 

El = bending stiffness, FL 

and 

n = 5 (pertinent variables) 
d = 3 (fundamental dimensions) 

N = 2 (size of complete set of dimensionless products). 


1) The n physical variables can be expressed in terms of 
the d basic dimensions (F, L, T) through the (d x n) 
matrix P 


F 

[Pi « L 
T 


m l 

1 0 
-2 1 
2 0 


AG El w 

1 1 0 

0 2 0 

0 0 -1 


The (n x d) matrix Z]_,and the (n x N) matrix can 
be arbitrarily selected to be 


1 

0 

0 


0 

0 

0 

1 

0 


0 

0 

0 

0 

1 

[z 1 = 

0 

0 

0 

0 

0 


1 

0 

0 

0 

0 


0 

1 


The (d x d) matrix P^ and the (d x N) matrix P^ are 


[P x ] 


[PI tz 1 i = 


’10 11 o " 


'1 0 0 “ 

= 

“10 1 

-2102 0 


0 10 


-2 10 

.2 0 0 0 -1 


0 0 1 

■ 

2 0 0 . 


0 0 0 



.0 0 0 . 



and 

[P 2 1 = [Pl [Z 2 1 = 


"l 0 1 1 o' 


0 o' 

sz 

"1 o' 

-21020 


0 0 


2 0 

.2 0 0.0-1 


0 0 


0-1 


1 0 



_0 1. 
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2) If the (N x N) matrix Q is defined to be 


IQI = 


1 o] , 

0 1 


the (d x N) matrix U is 


[U] = -[p l 1 [P 2 1 LQ1 


= 0 0 -f 1 0 


0 -1 -1 

-> • i 


[::] 


w 2 

-2 1 


3) The exponents are therefore 


[ V 


- [Z,l ID] + tz,l IQ] 


= 1 0 0 0 2 

0 1 0 -2 1 

0 0 1 -1 -j 

0 0 0 

0 0 0 


n 


C q ij J 


7T 1 n 2 
m f 0 l" 
2 

£ - 2 1 
AG -1 

El 1 0 

u 0 1 


4) The N = 2 dimensionless scaling parameters are 
determined from the vectors of q., as follows 


EI/AG£ 


m 1/2 £w/CAG) 1/2 or 7r 2 ' = m£ 2 io 2 /AG. 


M/8/7r//V MAM&TYA 

DENVER DIVISION 



3 -22 


MCR-68-87 


7T^ and 1^2 comprise the complete set of dimensionless 
scaling parameters and any other parameter can be 
obtained from them. For example 

* tt! A 2 
it * » _2_ = m& u) 

1 TT 1 El 
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4. VIBRATIONS OF BEAM-LIKE STRUCTURES 


Launch vehicles or certain portions of them, can be represented 
as beam structures. The modal characteristics are determined from 
the bending, torsional, shear, and longitudinal stiffnesses. Response 
is governed by the mass distributions, and the applied loads. Each 
of these quantities can be a function of position along the long- 
itudinal axis. 

Four classes of vibration problems are considered for dynamic 
scaling. The are identified as: (1) Free vibrations, no external 

loads; (2) Free vibrations, large external loads; (3) Forced 
vibrations, no external loads; (4) Forced vibrations large external 
loads. Notations, symbols, and arbitrary scale factors from Section 
3 are retained. For convenience we repeat the arbitrary scale factors 
and the resulting similitude requirements for force (F) , length (L) , 
and time (T). Subscripts m and p represent model and prototype 
respectively. 

n 

in 

n i = a 
p 

(la) 

E 

n„ = _tn 
E 

P 

(lb) 

^m 

n 0 * — 
3 P P 

(lc) 

After manipulating the dimensions of Z 9 E, and p we 
model requirements in terms of prototype dimensions from 

have for 
Equation 1. 

L - n- L 
m 1 P 

(2a) 

r 2 

F - n.. n„ F 

m 1 2 p 

(2b) 

T ■ n n ” 1/2 „ 1/2 T 

T m 1 2 3 T p 

(2c) 

Any of the physical properties of interest has one 
of these dimensions raised to some power. 

or more 

I i f' V, V R t) I V 1 5 1 O N 
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4.1 Free Vibrations. Mo External Loads 

It is assumed that beam rigidities are representative of 
cylindrical shells exhibiting beam-like behavior. End conditions 
for the model must be identical to those for the prototype. 

For the free vibration problem the dependent variable of 
interest is the frequency of vibration of the several mode shapes. 
Without external loads the frequency can be shown to depend on 
several independent factors such that 

W = f(pg, E, El, GJ, kAG) (3) 

Frequency w can be written nondimens ionally as 

7T » W& £. 

1 E 

From which we can write 


0 ) 


m 




w 

P 


(A) 


For this relationship to be valid Che following model terms 
must be scaled using preselected values of n^, n 2 and n^. 


P m g m “ n “ 1 n, p g 
mm x x p p 


EI - n. n_ E I 
mm 1 2 p p 


(5a) 

(5b) 


G J * n. n_ G J 
mm 1 2 p p 


kAG = n, n„ k A G 
m m m 1 2 p p p 


(5c) 

(5d) 


The most difficult modeling requirement to meet is given in 
Equation 5a. For B m “ g , both model and prototype in the same 
gravity field, the mass Hensity ratio requirement is shown in 
Fig. 4.1. The three lines show the effect of changing modulus. 

The upper line is a model material twice as stiff as the prototype, 
the middle line is the same modulus for model and prototype and 
the lower line is a model material 1/2 as stiff as the prototype. 

As the model becomes relatively smaller the model material must 
become increasingly more dense. 
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Fig, 4,1 Model Density Requirements as a 

Function Scale Factor for g = g 

m p 

By solving for n^ in Equation 5a and substituting in Equation 4 
the model frequency becomes 

co » n„ n_^^ 8 m oj for n. « n„n_-l ®m (6) 

m / 3 p. 1 l 3 *“ 

' ft 


If model and prototype are made of the same material and 
tested in the same g field, Equation 6 gives U) ■ w ; however, 
from Equation 4 for the same material m ^ 

-1 

u) ■ n, u . 
m 1 p 

To satisfy both requirements, it is apparent that n^ = 1 and the 
model becomes exactly the same as the prototype. 

For structures carrying no external loads or small external 
loads, one of the means of providing a different g environment 
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(for n. < 1 g > g ) is the use of a centrifuge. Model and 
prototype cfin ¥hen Ee modeled of the same material and the 
frequency relation, 

w ■ n, ^ w can be used 
m 1 p 


if 


-1 

g m * » g . 

m x p 


Another means of achieving the proper mass effect is to attach to 
the model additional material and not use centrifugal force. 

Added material must not provide rigidities that are not scaled 
or accounted for. 

When transverse bending modes are being studied, then Equation 
5b defines that rigidity requirement for the model. Torsional 
vibrations require the model condition given in Equation 5c. If 
large deformations are expected, or if the structure is short with 
respect to its depth. Equation 5d must be satisfied for the model 
to provide accurate dynamic response. 


Free lateral vibrations of beams can be found from the partial 
differential equations of the form 


El 



Pg 


a 2 
9 x 


3t y 


(7) 


Consideration of boundary conditions provides solutions for 
the response of a particular system. Exact geometric scaling 
provides the proper moment of inertia (I) and is accounted for 
by the scale factor n^. Material properties of the beam are 
included in the modulus term E, and scaled by the factor n^. 
Another property of the system is the unit mass density p and 
is included in the arbitrary scale factor n~. A term that is not 
arbitrary is the gravity constant g. The other two terms of 
Equation 7 are partial derivatives of the dynamic displacement x, 
fourth order in the length direction, second order in time, and 
are included in the modeling laws. 

4.2 Free Vibrations. Large External Loads 

For many of the practical cases the mass of the structure is 
small compared to other loads imposed on that structure. These 
may come from several sources such as payload, propellants, winds, 
movement through the atmosphere, control forces, trajectory 
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profiles, and others. The imposed moments, shears, torques and 
axial loads must be properly scaled for reliable model representa- 
tion of large scale launch vehicles. 

Continuing in the manner previously established, the frequency 
of a beam-like structure with a relatively small beam mass is a 
function of the following external loads and beam properties, 

0) - f(P, M, pg, Z, E, El, GJ, kAG) (8) 


The relationship for model and prototype frequency is 
-1 1/2 - 1/2 


“m ‘ "l n 2 


0) 


3 P 

Model requirements for the independent terms are 


(4) 


2 

P * n.. n„P 
m 1 2 p 

(9a) 

3 

M ■ n. n„M 
m 1 2 p 

(9b) 

p g » n ” n,p g 

mm 1 4 p p 

(5a) 

EI = n. 4 n„E I 
mm 1 2 p p 

(5b) 

G J - n. 4 n 0 G J 
mm 1 2 p p 

(5c) 

2 

kAG = n- n 0 k A G 
mmm 1 2 p p p 

(5d) 


The two vibration problems of one system with no external 
vibrating loads and another with large external loads is analogous 
to comparing a vibrating spring where in one case the spring only 
vibrates and the second case when the spring supports a large 
vibrating mass. In the first instance the inertial forces come 
from the vibrating spring and are important because they are the 
only inertial forces. In the second case when the mass is large 
compared to the spring mass the predominant inertial force is that 
of the mass. Scaling for this kind of load is accomplished using 
the similitude requirement of Equation 5a. Beam mass if sufficiently 
small is neglected and the supported mass oscillating in the same as 
the beam is scaled dynamically. 

Transverse vibrations are modified by axial loads on beams. 

A loaded launch vehicle on the stand or an accelerating vehicle 
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undergoes axial loading, and the magnitude of the load changes 
the frequency of vibration. Scaling this load is accomplished 
using Equation 9a. The model force is the square of the length 
scale factor times the modulus ratio times the prototype load. 

Another source of axial load, acting as a prestress, exists 
when launch vehicle tanks are pressurized. Lateral vibrational 
frequencies increase when longitudinal tension is applied to a 
beam and decrease when the beam is compressed. The tensile 
load from tensile prestress (pressure times cross sectional area) 
is scaled through the force term in Equation 9a. 

Shear loads, applied transversely to the axis of the beam, 
are scaled by the same factors used for axial loads. Equation 9a 
expresses the similitude requirements for shearing forces. 

Moments and torques are scaled by Equation 9b, the moment 
is applied in a plane containing the beam axis, and the torque 
in a plane normal to the axis. Whereas the force scale factor 
for the model is modified by the square of the length ratio the 
moment and/or torque on the model is a function of the length 
factor cubed times the modulus ratio. 

Lateral vibrations are resisted by the stiffness term El, 
scaled in Equation 5b. If this quantity is scaled and if the 
forces and masses that influence the response are properly scaled 
the free vibrations of the prototype structures can be found from 
the model from Equation 4. Torsional stiffness supplied by the 
GJ term provides the structural resistance to torsional vibrations, 
and the same frequency equation applies. 

When the model is an exact replica of the prototype and made 
of the same material, stiffnesses El, GJ, and kAG will be numerically 
correct and n~ and will be equal to unity. Force quantities 
will be scaled by n?, moment and torque by n.3 and inertial loads 
by n^ 1 , 

4.3 Forced Vibrations. No External Loads 

Beam structures not carrying masses, not subjected to axial 
loads, and being excited by time dependent forces can be dynamically 
modeled. Frequency becomes an independent variable, and is 
scaled according to the previously established similitude require- 
ments . Additionally, the maximum values of the exciting force 
applied to the model must be modeled assuming the prototype force 
is repetitive and of some known magnitude. 
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For the beam only the variables important to its response 
under forced vibrations with amplitude or the dependent variables 
can be given in the form 

A A 

x « f(w, P, g, £, E, El, GJ, kAG) (10) 

The hat superscript over u and P designate the frequency of the 
forcing function and the maximum value of the force. Any periodic, 
aperiodic, varying amplitude, or transient forcing system can be 
modeled if each of the frequency and load components of the complex 
system are scaled individually then combined and applied simultan- 
eously to the model. 

When complete similitude is achieved, the relationship between 
the response of the model and prototype is given by 

x ■ n., x (11) 

m 1 p 

Every component of the frequency spectrum of the forcing 
function imposed on the model is given by 

to - n- n. n, w , (12a) 

in x z j p 

and the corresponding force applied to the model when scaled 
dynamically becomes 

P - n/n, P n (12b) 

m 1 l p 

The p g term for unloaded structures is given by Equation 5a 
and applies ?o the forced vibration problem with the same discussion 
found for the free vibration problem with the same discussion 
found for the free vibrations. Remaining variables in Equation 9 
are again scaled as rigidity terms given by Equations 5a, 5b, and 
5c. 


4,4 Forced Vibrations. With External Loads 

External loads, not part of the forcing function, on a 
vibrating structure are of two kinds; constant with time, acting 
as a prestress; and inertial loads varying as mass times 
acceleration. Prestress loads exist as pressures, or dead loads, 
or acceleration loads and remain constant or vary slowly with 
time. Inertial loads can be varying because of changing velocity 
(acceleration) or changing mass (as a rocket) or both. 
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The terns for forced vibrations of beams carrying external 
loads are given by 

A A 

x « f(u), P, P, M, pg, J l t E, El, GJ, kAG) (13) 

Each of these terns has been previously scaled and no new similitude 
requirements are needed. However, for launch vehicles the external 
load P acting as a prestress either constant in time or varying time 
can be examined in greater detail. The definition of P comes from 
Newton's second law of motion so that 


P - ~ (pVv) 


To help establish the dimensionality of these terms, we take each of 
them individually to show using basic MLT dimensions, 

MLT~ 2 - “• (ML" 3 ) (L 3 ) (LT -1 ) 

and obviously after simple manipulation we have 

MLT -2 - MLT“ 2 

Therefore p is unit mass density of a given volume V, v is velocity, 
and the product of pVv is changing with time. For dynamic simi- 
litude from Equation 9a 

2 

P » n n 0 P 
m 1 2 p 


giving for the right hand side of Equation 14 


*77““ (p V v ) ■ n. n_ - —— (p V v ) 
dt ^m m m 1 2 dt P P P 

m p r v r 


The time rate of change of* momentum for the prototype must follow 
Equation 15. If the mass remains constant then 


(pV) - 0 
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dv 

v — 2L 

Pm m dt 

m 


n. 


n. 


P v 
P P 



(16) 


And if we designate 


dv 

dt 


as acceleration a, then from Equation 16 


2 Pp 2 -1 -3 

n l n 2 p V " n l n 2 n 3 n l 
m m 


or by combining terms the acceleration ratio becomes 


T “ n l n 2 n 3 (17) 

P 

This is an expected result when considering just the vibrating 
beams From Equation 5a 

-1 

p g » n. n 0 p g 
mm 1 2 K p°p 

Treating the g terms as acceleration terms and dividing by Pm 
the results of Equation 17 are readily found from Equation 5a 
above* 

If, however, the mass term is varying and the velocity remains 
constant, dv = q then Equation 15 becomes 
dt 


d 

V m dt 

m 


(P V ) 
m m 


2 d 

n i n n v TT" 
1 2 p dt 

P 


( Pn V r3 } 

P P 


(18) 


This describes the modeling for a system losing mass while 
maintining a constant velocity, and for launch vehicles this does 
not have much application* 


Most useful form of the equation for external forces on a 
rocket due to change of momentum is given by Equation 15 and the 
modeling laws can be applied to it in many forms* 


4 • 5 Loads From Fluids 


From the study of fluid mechanics we obtain a list of variables 
describing the physical quantities of gases* Not considering, at 
this time, temperature effects the several variables are listed in 
Table 4*1 using the engineering units and basic units. 
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Table 4.1 Quantities and Dimensions for Fluids 


Quantity 

Symbol 

Engineering 

Units 

Basic 

Units 

Force 

P 

F 

—2 

MLT 

Length 

z 

L 

L 

Velocity 

V 

lt" 1 

lt" 1 

Unit Mass Density 

p 

fl" 4 t 2 

ml" 3 

Dynamic Coefficient of 

14 

fl" 2 t 

ML-V 1 

viscosity 

Acceleration of Gravity 

g 

lt“ 2 

N> 

Speed of Sound 

c 

lt” 1 

lt" 1 

Surface Tension 

a* 

FL” 1 

CM 

1 

IS 


Nondimens ional numbers have been created by combining the 
quantities in Table 4.1, Similitude is achieved when the non- 
dimensional number has the same numerical value for the model 
and the prototype. See Ref, 3.6 for further discussion. 

For fluid simulation where viscosity is important the 
Reynolds' Number is scaled so that 

V Z p V i P 

m m m p p p 


Any combination of model velocity, length fluid density and 
kinematic viscosity that produces the same value as the prototype 
quantity satisfies this scaling requirement. 


To scale the resulting forces from fluid flow such as drag 
force the pressure coefficient relates force P to density, velocity 
and length so that 



o 2 0 2 n 2 0 2 

P v * P v * 

in m in p p p 


(19) 


Drag forces can be separated into two components, one 
viscosity (Reynolds' Number) and the other due to gravity. 


due to 
To 


*Not to be confused with stress in Table 3.1, 
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scale for gravity the Froude Number requires that 




( 20 ) 


When considering the effects of velocity in the sonic region or 
supersonic regime the Mach Number is given by 



( 21 ) 


For those fluids or scaling situations where surface tension 
effects are important then Weber’s Number defines a similitude 
requirement. To satisfy this requirement 


P v 
m m 


m 


m 


P v 

■ P.P, 


( 22 ) 


Although arbitrary scale factors such as £ / £ have not 
been defined for this section they can be applied En the same 
manner as for structures. Some of the forces P were not defined 
in the previous sections. In a launch vehicle these loads can 
come from many sources. When they occur as a result of actions 
between a structure and a fluid then considering the three basic 
dimensions of mass, length and time and by consistently modeling 
so that an arbitrary but unique relation exists 


L 

m 



L 

P 


M ■ 

m 




M 

P 


T - 

m 



T 

P 


and applying it to all physical quantities the model will give 
information regarding the prototype. 


4.6 Temperature Scaling 

When the temperature, or the quantity of heat, becomes a 
consideration for a prototype structure then other modeling laws 
and dimensions are necessary. The mechanical equivalent of heat 
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and energy establishes that a quantity of heat has the dimensions 
of work FL in engineering units and ML 2 !" 2 in basic units. 

Another dimension, temperature 6, and other definitions are 
introduced to express these physical quantities. These are 
given in Table 4.2 


Table 4.2 Quantities and Dimensions 
of Thermal Quantities 


Quantity 

Symbol 

Engineering 

Units 

Basic 

Units 

Temperature 

6 

e 

-l -l -l 

0 

-3 -1 

Heat Transfer Coefficient 

h 

k* 

FL T 0 A 

-1 -1 

mt e 

-3 -1 

Coefficient of Thermal 
Conduction 

ft e 

1 -9 -1 

MLT e 
2 -2 -1 

Specific Heat 

C 

H 

1 

CD 

LI 0 

-1 

Coefficient of Thermal 
Exp . 

3 

0 

0 


Several scaling numbers have become associated with heat 
transfer in gases. By using terms identified in Table 4.1 in 
combination with the thermal quantities in Table 4.2 we have 
from Ref. 3,6 


Grashof Vs 

Number 

« .-aas, a / 

y 2 

(23) 

Nusselt’s 

Number 

. in 

k 

(24) 

Praudtl's 

Number 

- & 

k 

(25) 

No Name 


- 

pv J ! 

(26) 


Another dimension, temperature 0, is introduced and is 
treated the same as length, mass, or time. The technique of 
scaling established for mechanical systems is simply extended from 
three dimensions to four dimensions. 


*k not to be confused with form factor k of kAG, 
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4.7 Nonlinear Effects of Structures 


Materials that are used to construct real structures and 
model structures can and do exhibit nonlinear behaviors. Some of 
the nonlinearities are exhibited by elastic and plastic behavior, 
strain-rate sensititivity, damping characteristics, and other 
phenomena observable on a macroscopic scale or in a microscopic 
snese. Generally when displacements are small compared to the 
least dimension of a structure, when frequencies are low, and 
when thermal effects are negligible, the nonlinearities are not 
important. 

For those cases where nonlinearities are important scaling 
is still possible if the nondimens ionalized response of the 
materials used for model and prototype can be plotted as the 
same curve. For example, when the nondimensional stress strain 
curves for model and prototype are the same curves, then the 
deflections of the model will scale as the length scale ratio, 
the forces will be changed by the length scale ratio squared, and 
the moments will be affected as the length scale cubed. 

The most important aspect of scaling is that a physical 
phenomena is not dependent on the numerical magnitude used to 
describe it. Physical laws when considered nondimensionally to 
relate model and prototype are the same physical laws and are not 
affected by size or dimension, all important physical quantities 
being considered. 
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5. SCALING OF SHELL STRUCTURES 


Shells used for launch vehicle structures may be monocoque 
construction of isotropic materials, orthotropic construction of 
isotropic materials, orthotropic construction of composite 
materials, or anisotropic construction of composite materials. 
Scaling for shell structures will be limited to isotropic materials, 
e.g. aluminum, titanium, magnesium and either monocoque (no 
stiffening) or orthotropic (stiffeners, rings) construction. 
Composite material or sandwich construction will not be considered 
although the scaling requirements are based on the same principles 
used for isotropic materials. 

Shell theory assumes the shell material to be an artificial 
substance with three moduli having infinite values and two 
Poisson ratios equal to zero. Using the coordinate system in 
Fig. 5,1 the several assumptions are given in Equation 1, 



Plane Tangent to Shell at 
Reference Surface 


Fig, 5.1 Coordinate System for Shell Material 

Using the directions identified in Fig, 5,1, the Hooke's 
law for shell materials takes the form 


; 11 


Vl + E 2 V 21 E 2 


( 1 ) 


°22 ' E l V 12 E l + E 2 E 2 


33 
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T 23 ' 0 


T 3i’° 


T = G Y 
12 12 '12 


This form of Hooke's law arises from the assumptions that 


E = G ■ G -*■ 00 
3 13 23 


( 2 ) 


V 13 " V 23 " ° 


e 3 * Y 31 " Y 23 


- 0 


These assumed values occur from the several hypothesis made 
regarding shell behavior such as normals remaining normal, no 
changes in shell wall thickness, small displacements, and thin 
shells . 

Although no material exists for which these assumptions 
are valid most materials used for shell construction approximate, 
sufficiently closely, this behavior. Model material for model 
shells made of materials other than used for the prototype may 
not be satisfactory if their properties are too different when 
compared to these assumptions. When approximating the shell as 
a beam with a hollow circular cross section, these considerations 
are unnecessary. When studying shell behavior these assumptions are 
required for tractable analytical solutions. Comparisons between 
experiment and theory are valid when both the test and the 
analysis use the same rules for material behavior. 

5,1 Monocoque Shells 

Dynamic scaling requires consideration of the three basic 
quantities of physical systems, mass, length, and time. These 
may be converted to the engineering units of force, length, and 
time. From previous developments (See Chapter 3) the three 
arbitrary scale factors will be chosen as 
l 

tn , _ . 

n . - T~ (3a) 
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5 ' 


n 


2 



n 3 



(3b) 


(3c) 


From Equation 3 the engineering units of length (L) , force (F) , 


and time (T) 

are modeled by using 


L 

» n.L 

(4a) 

m 

1 P 

F 

m 

" n l n 2 F p 

(4b) 

T 

ro 

1/2 -1/2 _ 

" n l n 2 n 3 T p 

(4c) 


Free vibrations of shells are associated with a circular 
frequency go, the dependent variable. Shell properties, material 
properties, and external loads are independent variables. Symbols 
and dimensions for these quantities are given in Table 3,1, A 
functional relationship can be written such that 

go - f (R, I, h* PSt E.» V, D, K, P, M, q) (5) 


A set of tt terms for these variables can be established as 
shown below. The dependent variable go is given by tt^ and the 
independent variables follow from on. 



For the independent factors from Equation 5 and recalling 
that Jt, E, and p have been chosen to define arbitrary scale 
factors given by Equations 3 the remaining terms are 

R 

n 2 " Z 
h 


(7a) 

(7b) 
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*5 

% 

"7 

*8 

TTg 

*10 


Pg 


1 


1* V. 


m 


pgR 

K 

2 

pgR 

p 

Ejl 2 


M 

" “ 

Er 


E 


(7c) 

(7d) 

(7e) 

(7f) 

(7g) 

(7h) 

(7i) 


Similitude requirements for Equations 7 to provide the 
equality of model and prototype it terms can be expressed as 


tt 

dl 

U 2p 

(8) 

7T q * 

3m 

% 


• 

« 

% 

• 

« 

• 

• 

• 


10m " 

*10p 



From the dimesnions of Table 3,1 and the relations expressed 
in Equations 4 the model dimensions become from Equations 7a and 
7b 


R 

m 


n- R 

1 p 


(9a) 


h 


m 


n, h 
1 P 


(9b) 


Model material properties and gravity requirements are found 
from Equations 7c and 7d 

0 
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®m - “l" 1 n 2 "S* 1 8 p * (9c) 

v m ' v p <9d) 


Bending rigidity D and extensional rigidity K are scaled 
as follows 


D 

m 


n. 


n D 

2 p 


(9e) 


K 


m 


“l n 2 K p 


(9f) 


External loads, moments, and pressures and internal loads, moments, 
and pressures are scaled for the model from corresponding prototype 
quantities by 

P m " n l 2 “2 P p (98 > 


M 

m 


n 2 M p 


(9h) 


\ " n 2 q p 


(9i) 


Equations 9a to 9i provide a working set of model requirements 
to find the dynamic response of a shell structure. If Equation 9b 
model thickness is scaled, then Equations 9e and 9f are satisfied 
when modeling with similar materials. If the structure Is unloaded 
or lightly loaded, its mass becomes important and (Eq. 9c) 
must be considered. Masses other than the structural mass that 
respond as the basic structure responds must also be scaled by 
making 


Applied loads, shears, moments, torques and pressures are 
made correct for the model by considering the modeling requirements 
imposed by Equations 9g, 9h and 9i. 


Nonlinearities known to exist in the prototype can be scaled 
into the model if the corresponding curves describing these non- 
linearities, when drawn nondimens ionally, are the same curve. An 
example of this is the elastic-plastic behavior shown by the 
stress-strain curve of a prototype material. If the nondimensional 
stress-strain curves for the modeling material and prototype 
material are the same curve then the model will respond in a scaled 
manner and will provide information about prototype response. 
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From 7J^ of Equation 6 the dependent variable of interest 
becomes for the model a frequency realted to the prototype by 


w 

m 





( 10 ) 


For models made of the same material as the full scale 
structure (n 2 * n^ ■ 1) the frequency relation is 



For dissimilar materials the terms n„ and n, for modulus and mass 
density affect the model response. Equation 10 applies only if 
the important modeling requirements outlined by Equation 9 are 
fulfilled. 

5.2 Stiffened Shells 

Many launch vehicles are made of stiffened shells with 
stringers running longitudinally. Combinations of rings and 
stringers are also used and in some cases rings alone may be 
used. Frequency response given by Equation 10 applies to the 
stiffened shell when all tt terms are scaled by Equation 9. For 
replica modeling the rigidity terms K and D are properly scaled. 

When replica modeling is not used, extensional and bending 
rigidities can be scaled. Expressions for these rigidities for 
stiffened shells are given in Table 5. When made nondimens ional 
by TTg and tt^ and made equal for model and prototype, then stiffnesses 
are scaled. Other variables must also be scaled as required. 

In Table 5.1 subscripts x and <j) refer to axial and circum- 
ferential directions respectively. Spacing of stiffening 
members is given by b^ measured in the axial direction and b^ 
measured circumferentially. Distance of the stiffener center 
of gravity from the reference surface (given by R) is given as 
dimension c. All other terms are as defined in Table 3.1, 
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Table 5.1 Stiffened Shell Rigidities (Continued) 
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6. SCALING OF LIQUID PROPELLANTS 


The dynamic behavior of launch vehicle liquid propellants 
is of considerable interest. At lift-off the propellants may 
comprise as much as 90% of the vehicle mass. The free surface 
of these propellants has a strong tendency to undergo large 
amplitude motion due to excitations caused by the flight environ- 
ment. Lateral accelerations arising from control surface deflections 
may contribute to these motions. 

With the development of large launch vehicles it has become 
difficult to provide adequate separation of slosh frequencies, control 
system frequencies and vehicle elastic frequencies. Knowledge of the 
dynamic slosh characteristics, including the coupled slosh-structural 
characteristics (e.g. , the well known POGO phenomonon) is of vital 
importance to loads and stability studies. An excellent treatment 
of the dynamic behavior of liquids in moving containers can be 
found in Reference 6.1. 

Scale models can be used to study lateral sloshing, coupled 
elastic-liquid oscillations, viscosity effects, compressibility effects, 
cavitational and surface tension effects. Typically, any model is 
designed for the simulation of some (but not all) of these phenomena. 

A general dimensional analysis of liquid sloshing phenomenon is 
developed (Section 6.1) for rigid and elastic tanks. Special limita- 
tions on several types of dynamic simulation of liquids are established 
(Section 6.2) and the results are presented in graphical form. These 
graphs may be used by the analyst to establish the overall simulation 
limitations if two or more liquid effects are considered simultaneously. 
The results of Section 6.2 are an extension of the results of Reference 
6 . 1 . 


6 .1 Simulation of Lateral Propellant Sloshing 

The investigation of the sloshing motions of liquid propellants 
in partially filled tanks is of importance to a complete dynamic analysis. 
If the slosh frequency is considerably below the fundamental structural 
frequency, slosh modes will couple with the flight control system and 
may produce instability. Slosh frequencies approximately equal to 
structural frequencies may cause considerable elastic-slosh coupling. 
Correct simulation of these coupled modes is essential to accurate 
loads and stability studies. 

Slosh modes can be suppressed to some extent through the use of ’ 
mechanical devices (e.g., baffles as discussed in Ref 6.2). It is 
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important to evaluate these devices through analysis or test, A 
mathematical approach is often not feasible due to the complexity of 
the boundary conditions introduced by the suppressing devices, however, 
the investigation of the problem through the use of suitably designed 
model experiments (Ref 6,3 and 6.4) is feasible. 

Simulation of slosh modes can be investigated for two general 
conditions: uncoupled and coupled with elastic modes. The uncoupled 

case can be considered when tank construction is such that the tank 
stiffness gives rise to structural frequencies which are not in close 
proximity to the lowest uncoupled slosh frequencies. On the other 
hand, there are many instances where the coupling of liquid modes and 
elastic tank modes becomes significant due to the proximity of the 
uncoupled natural frequencies. 

The important bending frequencies of current launch vehicles are 
usually less than 5 Hz and, as the important liquid surface modes in 
cylindrical tanks also occur at low frequency, coupling between the 
respective modes is likely. One of the few investigators to have con- 
sidered this phenomenon is Miles (Ref 6.5), This author developed the 
equations of motion for a circular cylindrical tank using a Lagrangian 
formulation and established the frequency equation in terms of uncoupled 
bending and sloshing frequencies and coupling coefficients. It was 
concluded that the liquid effects are only significant in those cases 
where the sloshing mass is a appreciable fraction of total mass. For 
typical launch vehicle structures the sloshing mass can be omitted 
from the calculation of bending frequencies and rigid- tank sloshing 
frequencies can be used in the determination of the tank bending 
characteristics. Results of an experimental investigation to confirm 
the analysis of Reference 6,5 are presented in Reference 6.6, 

A considerable amount of investigation has been devoted to the 
coupling of propellant motions with elastic tank breathing vibrations 
(e.g.. Ref 6.7 and 6.8). An extensive summary discussion of these 
efforts is presented in Reference 6.1. 

6.1.1 Uncoupled Propellant Sloshing 

The analysis of uncoupled lateral sloshing is based upon the 
following assumptions: 

a. Constant longitudinal acceleration, 

b. Sloshing arises from small amplitude excitations caused by 
either translational or rotational accelerations acting 
normal to the flight path. 
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c. Surface tension forces are negligible when compared to 
inertial and viscous forces, 

d. Forces due to gas pressure are negligible, 

e. Geometrical similarity between model and prototype. 
Translational Excitation 

A typical propellant tank subject to translational excitation 
is shown in Figure 6.1. 



Figure 6.1 Translational Excitation of Propellant Tank 

The resultant liquid force on the tank walls, expressed in 
functional form, is 

R = f (a, SL, X, y, T, p) 

where R = resultant liquid force on tank wall, F 

2 

a = longitudinal acceleration, L/T 
£ - characteristic length, L 
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A = any length (e.g., tank diameter (d), liquid depth (h), 
excitation amplitude (x q ) ) , L 

2 

y = liquid viscosity, FT/L 
p - liquid density, FT 2 /L 4 
t m excitation period, T. 

Application of the dimensional analysis techniques of Chapter 3 
yields four scaling parameters. These parameters are 

= RT 2 /pJl 4 , tt 2 = at 2 / SL, 

2 

tt 3 = p£ /yr and = A / £. 


The tt term is equivalent to the Froude Number and the tt term 
is equivalent to the Reynolds Number. Complete dimensional similarity 
is dictated by where the parameter A may be replaced by any 
desired physical dimension. 


If the analysis is developed using the resulting liquid pressure 
P as the dependent variable, the tt^ term becomes 

2 , 2 

TT^ = Pt /p£ 


which is a form of the Euler Number. The other non-dimensional terms 
are unchanged. 

Rotational Excitation 


The analysis is similar to that for translational excitation. 

The accelerations are determined by the location of the rotational 
axis and the amplitude of excitation. 

The functional relationship becomes 

M = f (a, £, A, u, p, T, 6 q ) 
where M » moment due to liquid pressure, F/L 

A = any length (including location of rotational axis (b ) ) , L 
0 q = excitation amplitude, 
and the other variables are as defined previously. 
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The five non-dimensional scaling relations are 

2 3 2 2 

= Mx / pJl , rr 2 = ax /£, tt 3 = p& /px» 

it = V A and tt c = 0 . 

4 DO 

As previously, tt^ is equivalent to the Froude Number and is 
equivalent to the Reynolds Number, Dimensional similarity is 
dictated by . Similarity of excitation amplitude follows from tt^ • 

The scaling relations for translational and rotational excita- 
tion can be examined for two conditions: 

a. Viscosity Neglected 

If viscous fluid forces are assumed to be negligible, the 
7T^ term becomes insignificant and the non-dimensional 
scaling relationship developed previously can be expressed 
as 


X /X 

= £ / £ = n, 

m p 

m p 

a /a 

= n(T /x ) 2 , 

m p 

p m 

R /R 

= n 3 (p /p ) (a /a ) 

m p 

m p m p 

P /P 

= n( -Pj a r) ’ 

r m *p m p 

m p 

H /M. 

= n 2 (p ,/pKa /a ) 
‘m ^p m p 

m p 

and 

CD 

3° 

li 

CD 

O 

no 


where the subscripts m and p refer to model and prototype, 
respectively. 

The first relation defines the geometric scale factor. The 
second indicates that any prototype acceleration can be simu- 
lated in a 1-g test environment through a time scale adjustment. 
The model- to-prototype frequency ratio, established by 7 r 2 f is 
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«/“ = t (a / a )/nl 1/2 . 

in p m p 


The expressions for force, pressure and moment indicate 
that the prototype fluid relations can be established 
following selection of the model fluid. 

b . Viscosity Included 

The scaling relations are identical to those presented 
above. In addition, the term must be satisfied. This 
relation can be written as 


n = (V m /v n )2/3(a n /a m )1/3 
m p p m 


where V - y/p is the liquid kinematic viscosity, L /T. 


This expression indicates that the length scale factor, n, 
is affected by both the density and viscosity of the model 
fluid. In order to investigate prototype accelerations at 
scale factors different from unity, the modeling liquid 
must have different properties than the prototype liquid. 

6.1.2 Elastically Coupled Propellant Sloshing 

The analysis of elastically coupled propellant sloshing is based 
upon the assumptions (a) through (e) of Section 6.1.1 but the complexity 
of the problem is magnified as the elastic structural tank properties 
must now be considered. The resultant liquid force acting on the 
tank walls, expressed as a function of the independent variables, is 


R = f (a, U V* v', \£, GJ, K x , B xy , B z> El, -f) 

where R = resultant liquid force on tank wall, F 

2 

a = longitudinal acceleration, L/T 
% = characteristic length, L 

X = any length (e.g. , tank diameter (d)* liquid depth (h), 
excitation amplitude (x. q ) ) 9 L 

2 

v = liquid kinematic viscosity, L/T 
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V 1 = Poisson’s ratio, 

= Poisson's ratio for bending, 

2 

GJ = tank torsional stiffness, FL 

K = tank ext'ensional stiffness, F/L 

= tank in-plane shear stiffness, F/L 

B = tank transverse shear stiffness, F/L 
z * 

2 

El = tank bending stiffness, FL 
T = excitation period, T, 

The non-dimensional scaling parameters, developed in accordance 
with the techniques established in Chapter 3, are 

TT 1 « Rv 2 t 2 /A 2 EI, tt 2 = at 2 / A, ir 3 = X/Jt» 

tt 4 = K x il 3 /El, ir 5 - vt/A 2 , tt 6 « v', 

*7 = v£, tt 8 - GJ/EI, 7T 9 = B z /K x „ 


and tt. n = B / K . 

10 xy x 

Note that tt ^ , TU and it, are similar to the scaling terns derived for 
uncoupled propellant sloshing. 

These scaling relations can be examined for two conditions: 


a. Viscosity Neglected 

If viscous fluid forces are assumed to be negligible the 
scaling relations which are dependent upon the length scale 
factor become 


X /X = 

m p 


l !% - n, 

m p 
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a / a 
m p 


= n(T / t ) 

P 


2 


and (K /El) /(K„/EI) = 1/n 3 . 

x m ^ p 


The first of these defines the model-to-prototype scale 
factor. The second indicates that any prototype acceleration 
can be simulated with an appropriate time scale. The 
frequency ratio, established by this expression is 

W /w » [(a /a J/n] 1 ^ 2 . 

m p m p 

Other scaling relationships, derived from TT,, 7Tg, it 

and tt are 
10 

C V EI) m /( V /EI) P ' ^ • 

(B z /EI) m /(B z /EI) p - 1/n 3 , 

(K /GJ) / (K /GJ) = 1/n 3 . 
x m x p 


The similarity requirements imposed by these other scaling 
equalities may be difficult to maintain when small scale 
models are fabricated. The inaccuracy of scaled skin 
thicknessess, due primarily to tolerance limitations, may 
cause significant deviations in these parameters , 


b. Viscosity Included 

Additional scaling relationships arise if viscous fluid 
forces are considered. The expression developed from uv is 

n - (v l/v ) 2 ^ 3 ( a t/ a m ^ 3 * 

m p pm 

which is identical to the scaling relationship of Section 
6,1.1, The expression for fluid forces, 

R /R = [(El) /(El) J/n 2 , 
tn p m p 

is established by and the above relation. 
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6.2 General Limitations on Liquid Scaling 

Density simulation or the scaling of slosh effects in a rigid, 
uncoupled tank (Section 6.1.1) presents no difficulties. Scaling 
of liquid propellants to other conditions may impose limitations on 
physically feasible models. Typical scaling requirements include 
gravitational-viscosity scaling, gravitational-compressibility 
scaling, gravitational-cavitation scaling and gravitational-surface 
tension scaling. 

The non-dimensional scaling relations for each of the above 
cases can be established through examination of the functional 
relationship 


f(AP, £, a, v* c, p, r, c/g, 0) = 0 

where AP - P - P (P = gas pressure above liquid, 

gas vapor gas 2 

^vapor = va P or P ressure of liquid), F/L 
£ = characteristic length, L 

2 

a = longitudinal tank acceleration, L/T 

2 

V = liquid kinematic viscosity, L /T 

c = velocity of sound in liquid, L/T 

2 4 

p = liquid density, FT /L 

r = radius of meniscus in equivalent capillary tube, L 

a = liquid gas surface tension, F/L 

2 

g = acceleration of gravity, L/T 

0 = contact angle between liquid and tank wall. 

It is assumed that the solid-liquid and solid-gas interface effects 
are negligible. 


Application of the dimensional analysis techniques of Chapter 3 
yields six scaling parameters* These parameters are 

tt ]L - £ 3 a/v 2 , tt 2 = c 2 /a£, tt^ = AP/pa£, 

2 

tt^ = r /£, TT^ = pgr /o and n = 6 . 
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The scaling requirements imposed by these parameters and the 
limitations on scaling of liquids in a rigid tank are discussed in 
Sections 6.2.1 through 6.2.5. Gravitational-density scaling of 
liquids in elastic tanks is discussed in Section 6.2.6. 


6.2.1 Gray i tat io nal~ V is co s i ty Scaling 


The significant scaling parameter for gravitational-viscosity 
simulation is 

= 9? a/V 2 . 


A bar chart which compares the kinematic viscosities of typical 
propellents and possible model liquids is presented in Figure 6.2. 
The kinematic viscosities are shown for the temperature ranges which 
may be experienced during launch vehicle boost flight. 



0.002 0.005 0.01 0.02 


Figure 6.2 Kinematic Viscosity, Typical Missile Propellant and Model Liquids 
(Based on Fig. 5.4 and 5.5, Ref 6.1) 
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and the subscripts m and p refer to model and prototype, respectively. 

The relation between the these parameters is presented in 
Figure 6.3. 



Figure 6.3 Simulation Chart: Gravitational-Viscosity Scaling 

(Based on Fig. 5.6, Ref 6.1) 
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The difficulty which arises if gravitational-viscosity simulation 
of cryogenic propellants is attempted becomes evident. The kinematic 
viscosity ratio must be approximately equal to or greater than unity 
which makes it virtually Impossible to preserve the required scaling 
parameters, for typical boost acceleration ranges, with model scale 
factors less than unity. 

6.2.2 Gravitational-Compressibility Scaling 

The controlling scaling parameter for gravitational-compressibility 
scaling is 

2, A 

^2 ■ c /aJl. 


The required model-to-prototype scaling condition, established 


by i\ , is 


c 2 /a l - c 2 /a SL 
m mm p p p 


or 


c = (na ) 
r r 


1/2 


where c « c /c . 

r m p 

The relation between these three parameters is presented in 
Figure 6.4. As typical values of the sonic velocity ratio will lie 
in the range 1/4 < c < 4, gravitational-compressibility propellant 

scaling for typical Soost acceleration ranges is possible for a 
limited range of desirable model scale factors. The simulation 
of nitrogen tetroxide or LOX is difficult (at desirable scale 
factors) as the sonic velocity ratio, assuming the use of a common 
modeling liquid, will not be much less than unity. 

6.2. 3 Gravitational-Cavitation Scaling 

The scaling parameter for gravitational-cavitation scaling is 

tt 3 = Ap/p a l. 

This scaling requirement can be satisfied, for any length 
scale and prototype liquid, by adjusting the gas pressure (assuming 
that the absolute gas pressure is not determined by any other scaling 
parameter). The ratio AP/p, for typical model liquids and an assumed 
model gas pressure of one atmosphere, will lie between 0.6 and 1.4. 

The required model-to-prototype scaling condition, established by ir , is 
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(AP/p) /aA = (AP/p) /a. l 
m m m P P P 

or (AP/p) r = na r 

where (AP/p) = (AP/p) / (AP/p) . 

r m p 



(Based on Fig, 5.7, Ref 6.1) 
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The relation between these three parameters is presented in 
Figure 6.5. 



Simulation (at atmospheric pressure) of a cryogenic prototype 
fluid or a fluid with a very high vapor pressure does not appear to 
be feasible for convenient model scale factors. Boost acceleration 
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simulation will require model gas pressures to be reduced below 
atmospheric. Gravitational-cavitation effects in a low-g 
environment can be simulated with small scale models and 
atmospheric gas pressures. 

6.2.4 Gravitational-Surface Tension Scaling 

Gravitational-surface tension scaling is controlled rr^ , 
and Tig of Section 6.2. tt^, a variation of the Bond Number, can 
be combined with tt^ to yield 

tt 5 * - pg H 2 /a. 


The required model-to-prototype scaling, expressed in terms of the 
kinematic surface tension, (J) = c/p, is. 


or 


g & /4> 
m m 




2 

n a 


r 


where <f> = <J> / 4> and 

r m p 




a /a 
m p 


a a 


r 


The relation between these parameters is presented in Figure 6.6. 
Velues of <j) , for typical prototype propellants (e.g . , kerosene, LOX) 
and model fluids will lie between 0.1 and 10. It is obvious that 
surface tension scaling with desirable size models at small accelera- 
tion ratios is not feasible. The simulation of low gravity surface 
tension effects is possible with small scale models in a 1-g test 
environment . 

The scaling relation imposed by tt^. 



is probably insignificant for dynamic simulation as the contact angle 
influences interface area variations in the energy equations as cos 0. 
Nevertheless, this requirement eliminates the possibility of simulating 
a wetting fluid with a non-wetting fluid and vice versa. 

The discussion of Sections 6.2.1 through 6.2.4 has been restricted 
to the simulation of inertial effects with any of four of the more impor- 
tant liquid effects. The results of these sections may be used concurrently 
if the analyst wishes to consider the simultaneous simulation of several 
of these effects. 
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Figure 6.6 Simulation Chart: Gravitational-Surface Tension Scaling 

(Based on Fig. 5.10, Ref 6.1) 


6.2.5 Relative Importance of Capillary and Gravitational Effects 

The dynamic characteristics of a gas-liquid system may be 
influenced by several forces (e.g., viscous forces, capillary forces, 
body forces) but often all but one of these forces can be regarded as 
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negligible. The hydrodynamic behavior of the system can be 
separated into approximate regimes through examination of several 
dimensionless parameters. 

2 

The Weber Number, We « pV £/cr, defines an estimate of the ratio 
of inertial to capillary forces. This parameter appears (in a 
slightly different form) in Section 6.2.4 as the non-dimensional 
scaling relation tt,.*. Note that for We much greater than unity the 
capillary forces are insignificant. For We much less than unity the 
opposite is true. 

2 

The Froude Number, Fr ■ V /g£, relates the inertial forces to 
the body forces. It is evident that body forces play an important 
role in the fluid dynamics for values of Fr much less than unity. 

The regimes defined by the Weber Number, Froude Number and Bond 
Number (Bo = We/Fr) are indicated in Figure 6.7 taken from 
Reference 6.1. 



Figure 6,7 Hydrodynamic Regimes (See Fig. 11.5, Ref 6.1) 


This figure gives no indication of the relative importance of viscous 
effects and the effect of these forces must be established separately* 
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An order-of-magnitude estimate of the time required for liquid 
reorientation following the transition between two hydrostatic regimes 
can be obtained through examination of the "characteristic response 
time." An analysis which develops expressions for these times is 
presented in Section 11.4 of Reference 6.1. The characteristics times 
are determined to be 

1/2 

t = (5,/g) for the gravity dominated regime 

3 2/2 

and t * (p£ /C) for the capillary dominated regime. 

Figure 6.8 indicates characteristic response times for the various 
hydrostatic regimes. 


1000 
100 
10 
1 

0.1 

10 ' 8 1 0" 6 1 o " 4 10" 2 1 100 

Acceleration, g 

Figure 6.8 Characteristic Response Times (Based on Fig. 11.6, Ref 6.1) 
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6.2.6 Elas tic-Gravitationai-Densi ty Scaling 

The simulation of an elastic launch vehicle structure containing 
a sloshing liquid has been developed in Chapter 4 where it was assumed 
that significant coupling exists between the slosh modes and the beam- 
like structural modes. The controlling scaling parameters, established 
in Chapter 4, are 

= m Aj 2 /EI, tt 2 - Sm 2 / a and TT = p$, 2 /m 

2 2 

where m = mass/unit length, FT /L 

H = characteristic length, L 
co = frequency, 1/T 

2 

El = beam stiffness, FL 

. 2 

a - acceleration, L/T 
p = liquid density, FT /L* 

and geometrical structural similarity is required. 

The acceleration ratio, using the subscript r to denote model- 
to-prototype ratios, is 

a = (El) /m n 3 
r r r 

where n = & / !L t The mass ratio, due to the assumption that the 

ratio of sloshing ?o non-sloshing masses be identical for model and 
prototype, is 



and the acceleration ratio then becomes 

a = (El) / P n 5 
r r r 

Replica scaling requires geometric similarity and idential 
materials. This implies that ■» 1, I = n^ and P * 1 so that the 
acceleration ratio becomes 
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a 


r 


*» 1/n 


which clearly indicates the impossibility of modeling the boost 
acceleration phase with small scale replica models. However, the 
possibility of using dissimilar materials must not be precluded 
although tank geometric similarity is always required to preserve 
adequate slosh simulation. If it is assumed that 


i 3 t 
r r 


n 3 t , 
x 9 


where t defines the ratio of model-to-prototype tank wall thickness, 
the acceleration ratio becomes 


a 


r 


E t /p n 
r r r 


2 


A simulation chart in terms of the parameter Q = E t /P is 
presented in Figure 6.9* 

c 

<u 
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As typical model sizes and boost phase accelerations yield 
values of Q ranging from approximately 0.03 to 0.002* the 
applicability of metal models is limited. The liquid density ratio 
(unless mercury is used as the model liquid) cannot exceed about 
2.0; the elastic modulus for aluminum prototype and magnesium model 
would be 0.6 and a thickness ratio of 0.1 would be required to 
achieve Q - 0.03. This ratio of tank wall thicknesses may create 
seriois manufacturing* assembly and handling problems. To achieve 
lower values of Q requires the use of plastics in model construction 
and/or the use of mercury as the model liquid. 

6.3 Simulation of Non-Linear Longitudinal Effects 

The effect of longitudinal excitations upon liquids depends to 
a great extent upon the amplitude and the frequency of the excita- 
tions. At relatively low frequencies the free surface of the liquid 
may respond in a large amplitude standing wave with a frequency 
equal to one-half the excitation frequency. At high frequencies the 
liquid amplitude is diminished but the waves may disintegrate and 
generate a low frequency standing wave of large amplitude. At very 
high vibration input levels small vapor bubbles can become 
entrapped in the liquid. These bubbles can achieve negative 
buoyancy and sink to the tank bottom. 


A historical survey and a summary of analyses and conclusions 
pertaining to these non-linear effects has been presented by 
Abramson (Ref 6 . 1) . 


6.3.1 Subharmonic, Superharmonic and Other Non-Linear Excitation 

Free surface standing waves are generated when a cylindrical 
tank containing a liquid is forced to oscillate vertically. The 
primary liquid response occurs at one-half the excitation frequency 
and the liquid motion is usually identified as the one-half sub- 
harmonic response. Other response modes may occur at the excitation 
frequency as well as at superharmonic frequencies. This phenomenon 
is essentially non-linear and is a suitable area for simulation 
studies . 

The functional relationship between liquid pressure and other 
independent variables in an assumed rigid tank is 

p ■* f(g, K p» w, fc) 

2 

where F - liquid pressure at depth d* F/L 

2 

g = steady state tank acceleration, L/T 
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X •> any length (e.g., excitation amplitude (x ), 
radius of tank (R), liquid depth (d)), L 

P * liquid density, FT^/L^ 

w =» excitation frequency, 1/T 

2 * characteristic length, L. 

Application of the principles of dimensional analysis yields 
three non-dimensional scaling relations. These relations are 

Ttj^ » 0Q 2 Vg, ^2 = and ^3 ■ x/st,. 

The first term, when written as 
TT^ • to x Q /g 

defines the dynamic-to-static acceleration ratio, may be 
expressed as 


"*2* " p /Pg d 

to define the dynamic-to-static pressure ratio. Complete geometric 
similarity is dictated by 

6.3.2 Bubble Cluster Formation 


Another non-linear phenomenon associated with vertical tank 
excitations is the formation of bubble clusters. These clusters arise 
due to the apparent negative buoyancy of small gas bubbles which, 
above some limiting oscillatory acceleration level, become trapped 
in the liquid. This phenomenon results in the periodic growth and 
destruction of gas bubble clusters below the liquid surface. 

The functional relationship between liquid pressure and other 
variables is 


P = f (g,. X, p, w, SL, c, P q ) 

where c - sonic velocity in liquid-bubble mixture, L/T 

2 

» gas pressure, F/L 

and the other variables are as defined in Section 6.3.1. 


C0/?P0l?/ir/0/V 

DENVER DIVISION 



MCR-.68-.87 


6-23 


The five scaling relations developed from the above are 
ir ]L = wA/c, tt 2 = X/ £, tt 3 = P Q /pgil, 

2 

= a) X/g and = P/pgJi. 

when expressed as 
TT^* = cod/c, 

defines a dimensionless pressure wave length, TT- can be written 
as 

^3* = P c/ P8d 

to define the gas-to-liquid static pressure ratio and tt 5 as 
tt 5 * = P/pgd 

to define liquid dynamic- to-static pressure ratio. The dynamic-to- 
static acceleration ratio follows from 

V = “V 8 

and geometric similarity is dictated by tt 2 . The ratio of gas to 
liquid volumes 


*6 ■ Y v * 


follows from T^. 
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7 . DISTORTED SCALING OF LAUNCH VEHICLE STRUCTURES 


Stiffened cylindrical and conical shells are typical of tank 
walls, interstage adapters and fairings. Proper simulation of their 
dynamic characteristics is required if accurate simulation of the 
overall launch vehicle dynamic characteristics is to be realized. 
Frequently, the thin skins that form a major portion of these 
structures preclude the use of replica scaling (Chapter 9). Replace- 
ment of the prototype shell structure with an equivalent distorted 
model is then necessary. 

There are several approaches to the analysis of shell configura- 
tions. Classical solutions are available for most axisymmetric shell 
geometries (e.g.. Ref 7.1) but many practical problems involve 
unsymmetrical or stiffened shells for which classical solutions are not 
readily available. Due to the complexity of the mathematical formula- 
tion, stiffened cylinders and cones are generally investigated through 
the use of finite element techniques. These techniques idealize the 
actual shell and stiffeners as structural systems of plates and beams 
having equivalent stiffness and mass. Application of these methods 
to symmetrically and unsymmetrically stiffened cylinders and to 
symmetrically stiffened cones is presented in Reference 7.2. 

A somewhat analogous method of distorted shell modeling which 
appears to be applicable to typical launch vehicle structures is 
truss-ring modeling. This technique differs from that outlined in 
Reference 7.2 in that the shell is replaced by truss-ring modules 
consisting of transverse frames connected by a framework of axially 
loaded bars. Application of these methods to the dynamic modeling 
of shells is discussed in Section 7.1. Simulation of orthotropic 
cylinders using reinforced plastic models is discussed in Section 7.2. 

7 .1 Truss-Ring Modeling of Cylindrical Shells 

The cylindrical shell considered here is a continuous three- 
dimensional skin structure. It is proposed that the shell be 
suitably divided into a number of units, each unit consisting of a 
number of segments. Therefore, the entire cylindrical shell can be 
considered as an assemblage of a finite number of segments. 

The extremely thin skin gages which may result from replica 
scaling dictate that the cylindrical units be simulated with a distorted 
model. The model suggested to represent the continuous elastic 
properties of a typical cylindrical unit is defined as a "module." 

This module (Fig, 7,1) consists of an elastic ring and several 
axially stressed bars. The bars are equally spaced about the ring 
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circumference and are singly connected at each end by frictionless 
pin joints. These joints attach to the module ring on one side and 
to the ring of the adjacent module on the other side. Note that 
the geometrical arrangement of the truss bars must be symmetrical 
about any ring diameter. 



Truss Bar 


Figure 7,1 Truss-Ring Module 


The proposed module must satisfy certain preassigned static and 
dynamic conditions. The elastic ring will simulate the radial 
effects while the bars are assumed to take only axial forces. The 
combination of truss and connecting ring provides the required overall 
simulation of longitudinal, lateral and torsional stiffnesses. The 
entire cylindrical shell may be simulated by a number of these modules. 
As the bars are pin- connected, the truss between adjacent rings forms a 
discrete system. Thus, the assemblage of the several modules is 
precisely an assemblage of discrete units and the entire truss-ring 
modeling analysis can be deemed to be analogous to a finite element 
approach . 

The following sections detail the equivalent longitudinal, 
torsional, bending and lateral stiffnesses of a module. Buckling 
simulation is not considered in this analysis. The simulation of 
shell buckling characteristics is a logical extension to these 
analyses . 
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7.1.1 Geometry of Module Components 

A typical module (Fig. 7.2) consists of an elastic ring and a 
number of pin-connected bars. 



Figure 7.2 Geometry of Module Element 


The number of joints attached to a ring is denoted by n and therefore 
the number of bars is equal to 2n, Geometrical symmetry dictates 
that all bars have identical length. 

The notation used in the subsequent analyses of module stiffnesses 

is: 

R = radius of ring centerline 

n = number of joints on one side of ring, 

d = longitudinal distance between ring centerlines. 
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2 2 2 

£ = length of truss bar, £ * a + d , 

a « n/n, 

0 b cos ^(d/£), 
a *= 2R sin (ct/2). 

A set of x, y, z axes are employed such that the x-axis defines 
the longitudinal module direction and the y and z axes define the 
lateral directions (Fig. 7.2). The axial bar force deviates from the 
module longitudinal axis by the angle 0 and may therefore be decomposed 
into x, y and z components for the individual stiffness analyses. 

7.1.2 Longitudinal Stiffness of a Module 

The longitudinal module stiffness is defined as the longitudinal 
force required to produce unit longitudinal module deformation. The 
connecting rings are assumed rigid. 

The length of a typical rod, expressed in terms of the module 
length d and chord distance between joints a, is 

£ 2 - d 2 + a 2 . (1) 

The relationship between the variation in rod length <5£ with 
change in module length <Sd can be obtained from the differentiation 
of Equation 1 and is 


6£ = (d/£) <5d (2) 

The assumption of a rigid ring dictates that the dimension a remain 
constant. Assuming that the deformation is within the elastic 
range so that Hooke's Law is applicable allows the axial rod force to 
be expressed as 

F i = k(<S£) - k(d/£) 6d (3) 

where the axial stiffness of an individual rod is 

k = AE/£. (4) 

It is to be noted that no subscript is used in Equation 4 as the 
rod area A, length £ and material elasticity E are assumed identical 
for all rods. 
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The longitudinal component F^ of the axial rod force can be 
computed from Equation 3 and is 

F ix = F i (d/Jt) = k W/^) 2 (5) 

Each module consists of 2n rods and the total longitudinal module 
force is 

\ = 2nF = 2nk(d/&) 2 fid. (6) 

The longitudinal module stiffness, expressed in terms of the rod 
stiffness and the module geometry, is 

k L - F /fid - 2nk(d/&) 2 . (7) 

In model analysis it is convenient to employ the concept of 
an equivalent elastic parameter, For longitudinal extension the 
equivalent longitudinal elasticity AE of a module can be defined 
as 

AE = kjd = 2nAE(d/£) 3 . (8) 

Note that no forces but the longitudinal resultant force of 
Equation 6 are required for unit longitudinal module deformation. 
This is a consequence of the assumptions of rigid ring and module 
symmetry. As the applied force and resulting deformation act in 
the same direction, the longitudinal stiffness is defined as a 
direct stiffness. The lateral components of the axial rod forces 
are mutually cancelled due to the symmetrical connections of truss 
loads to ring joints and no cross— stiffness terms exist. 


7 .1.3 Torsional Stiffness of a Module 

The torsional module stiffness is provided by the axial tension 
and compression of the truss rods. Again it is assumed that the 
ring is rigid so that only a torque is required to produce twisting 
deformation. The torsional stiffness is defined as the torque 
required for unit twist angle and can he derived by rotating the 
rigid ring about the longitudinal axis through the angle 60 . The 
pin-joint attached at the ring circumference will travel an arc 
distance 6a (Fig. 7.3) where 

6a = R(69). (9) 
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Figure 7 .3 Torsional Displacement of Truss Triangle 


Since the change in rod length SSL is a function of 6a and the module 
length d is invar lent, the differentiation of Equation 1 yields 

6 A- R(a /SL) 69. (10) 

The axial rod force F^, assuming that Hooke's Law is applicable, is 

¥ ± = k(6A) = RR(a/ SL) 66 (11) 

where k is as defined in Equation 4. This force can be decomposed 
into two components; one in the longitudinal direction and the other 
in the plane of the ring. The latter, which produces a twisting 
moment about the center of the ring (Fig. 7.3) can be expressed as 

F ia “ * kR(a/A) 2 66. (12) 
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The twisting moment developed by this force is 

T = F. R cos (a/2) = k(Ra/Jt) 2 cos (a/2)S6. (13) 

i ia 

There are two rods associated with each pin joint; one in 
compression, the other in tension. The longitudinal components of 
the two rod forces are mutually cancelled but the in-plane components 
F^ a are equal in magnitude and act in the same projection direction. 
Hence, the total twisting moment T is 

T = 2n T ± = 2nk(aR/&) 2 cos (a/ 2) 50. (14) 

and the torsional stiffness is 

kj, = T/60 = 2nk(aR/A) 2 cos (a/2). (15) 

The equivalent torsional elasticity is 

GJ = k T d = 2jiAE(d/i) (aR/£) 2 cos (a/2). (16) 

The cancellation of the longitudinal rod force components at 
an individual joint implies that the torsional stiffness is 
independent of the longitudinal stiffness. Due to the assumed 
symmetry of the rods, the twisting force components provide neither 
sectional lateral shear or bending deformations. The torsional 
stiffness is therefore a direct stiffness. 

7.1.4 Bending Stiffness of a Module 

The bending stiffness is defined as the bending moment required 
to produce a unit ring rotation about a diametric axis (Fig, 7.4a). 
The ring is assumed to be rigid. The influence of rod forces will be 
discussed in a later section. 

The truss rods are assumed to be arranged symmetrically about 
any diameter and it follows that two neutral axis positions must be 
considered. The bending neutral axis may either pass through the 
joints (Fig. 7.4b) or lie between the joints (Fig. 7.4c). 

In order to establish the moment required to produce 6<|> it is 
necessary to consider the axial deformation of the truss rods. Let 
the joints be numbered counterclockwise as shown in Figure 7. 4b, 

The initial joint is defined by the angle a and the angle between 
adjacent joints is 2a, The jth joint is then defined by the angle 

= 2ja - a o (17) 
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where a » 0 if the neutral axis is as shown in Figure 7.4b and 
a = a if the neutral axis is as shown in Figure 7.4c. The 
d?stance of the jth joint from the neutral axis is 

- R sin 3^ - R sin (2ja - a Q > (18) 

and the longitudinal displacement x^ is 

X j * y j C 19 > 

It follows that the axial rod deformation at the jth joint is 

« yj(dM)<5<j> (20) 

and the axial rod force, assuming a linear elastic material, is 

Fj *> k = ky^(d/£) 6cJ> (21) 

where k is as defined in Equation 4. The longitudinal component 
of Fj la 

F - FjCd/l) - ky^d/l) 2 6<|> (22) 

and, as there are two identical components at the jth joint, 
the total longitudinal joint force is 

2F jx = 2k yj (dM) 2 6c}>. (23) 



Figure 7.4 Bending Displacement of Truss-Ring Module 
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The bending moment about the neutral axis produced by the total 
joint force can be expressed as 

= 2*^^ = 2k(yjd/A) 2 6<() (24) 

and the total moment is 


n 

n/2 

r n/2 .. 



-E 

M 3 ■ 2 L “3 

« 4k(d/n) 2 y/ 

6<f> 


j=i 

j-1 

L j«=i 





- n/2 




- 4k (Rd/ Jl) 2 

^ ^ sin 2 (2ja - a ) 

6<j>a 

(25) 



- j»l J 




The sum may be taken over half the circumference because the moments 
produced by the joint forces are symmetric about the neutral axis. 

The bending stiffness is 

n/2 

kg = M/6<t> = 4k(Rd/A) 2 ^ sin 2 (2ja - a Q ) (26) 

j=l 

and the equivalent bending elasticity can be expressed as 

n/2 

El = k D d = 4AE (d/ Z) (Rd/Z) 2 > sin 2 (2ja - a) . (27) 

IJ D O 

j=l 

The bending stiffness is a direct stiffness. However, it 
should be noted that this bending stiffness couples with a lateral 
cross-stiffness as a lateral force is required to prevent lateral 
ring displacement under the influence of the bending moment. 

The formulation of the lateral cross-stiffness is omitted. 

7.1.5 Lateral Stiffness of a Module 

The lateral stiffness is defined as the lateral ring force required 
to produce unit lateral displacement. All other displacements are zero 
and a rigid ring is assumed. 
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Consider the half circular ring and connecting joints shown in 
Figure 7.5. 



Figure 7.5 Lateral Displacement of Truss-Ring Module 


The ith joint attaches to the ring and the jth joint attaches to 
the adjacent ring. It is evident that the elongation of the ith 
bar due to lateral displacement, Sy, is 

6a ^ = 6y cos (8^ - &/2) (28) 

and the compression of the jth bar is 

6a ^ = <5y cos (8^ - a/2). (29) 

If the truss rods and joints are numbered as shown in Figure 7.5, 
either Equation 28 or Equation 29 may be used to express individual 
rod deformation provided that 8 is restricted to 

0 < 8 < tt/2 (30) 
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The undeformed rod length has been given previously (Equation 1). 
As it is assumed that the module length d is held fixed, differentia- 
tion of Equation 1 yields 

2 26£^ = 2a 6a^ or <$£^ * (a/ £) cos (8^ - a/2)6y. (31) 

The axial bar force, assuming that the deformation is within the 
elastic range so that Hooke's Law is applicable, is 

=» k 62^ » k(a/ 2) cos (8^ “ a/2) 6y (32) 

where k is as defined in Equation 4. This force can be decomposed 
into two components; one in the longitudinal direction and the other 
along the projection direction (Fig. 7.5). This latter component is 

F ia = F^(a/ 2) » k(a/&) 2 cos (8 ± - a/2)6y. (33) 

The lateral component of F is then 

■F - F^ cos (8 a - a/2) ■ k(a/2) 2 cos 2 (8 A - a/2)<Sy (34) 

and the total lateral force, determined by summing the rod forces 
over the quarter-circle defined by Equation 30, is 

n/2 -n/2 

F y = 4 ^ F iy = 4k (a/ 2) 2 ^ c.os 2 ('$ 1 - a/2) 6y. (35) 

i=l L i«l 

The lateral stiffness can now be expressed as 

n/2 

k y = F / % = 4k (a/ A) 2 ^2 cos 2 (8^ - a/2) (36) 

i=l 

and the equivalent lateral elasticity as 

n/2 

EI l = kd - 4AE(d/2)(a/2) 2 ^ cos 2 (0 ± - a/2) . (37) 

i=l 

The lateral stiffness is a direct stiffness. However, it should 
be noted that this lateral stiffness couples with a bending cross- 
stiffness as a bending moment is required to prevent rotational 
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ring displacement under the influence of the lateral force. The 
formulation of the bending cross-stiffness is omitted. 

7.1.6 Coupling Between Circumferential and Longitudinal 
Deformation 


The analysis of module longitudinal stiffness (Section 7.1.2) 
is based upon the assumption that the truss rings are rigid. If 
the ring is considered to be elastic the effect of the axial rod 
force components which lie in the plane of a truss ring must be 
considered. 

Consider a module subjected to a longitudinal force F. The 
longitudinal force component in an individual rod is 

F ix = F/2n (38) 

and the component along a chord projection in the ring plane is 

F. a - (F/2n) (a/d). (39) 

This force component can be further decomposed into two components * 
One of these components is tangential to the ring; the other is in 
the radial direction. The tangential components of two rods which 
connect at a single joint are mutually cancelled but the radial 
components introduce circumferential deformation (Figure 7.6a). 

The radial force component in a individual rod is 

F ir = F ia Sin (0i/2) = (F/2n).(a/d) sin (a/2) (40) 

and it follows that the total radial force acting on the elastic 
ring can be expressed as 

F r = 2nF ir = F (a/d') sin (a/2) . (41) 

This resultant force can be approximated as a distributed 
compressive force acting about the circumference of the ring if 
the number of truss rods is large or the chord projection length 
is small. This approximate uniform force per unit length about the 
circumference is equal to 

P r = F r /2TTR * (F/2 ttR) (a/d) sin (a/2). (42) 
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(a) 

Figure 7.6 Circumferential-Longitudinal Coupling 


For a uniform thin ring the circumferential force F (Fig. 7.6b) is 
obviously equal to c 

F c = P r R = (F/2TT) (a/d) sin (a/2) (43) 

and the unit circumferential strain in the ring can be expressed as 

£ c = V A r E r = < F / 27TA r E r > (a/d) sin (a/2) (44) 

where and define the ring elasticity and cross-sectional area. 
This strain should be equivalent to the circumferential strain of the 
simulated shell. 

The change in radius, AR, is equal to the total circumferential 
deformation divided by 2 tt, or 

AR = 2ttR£ / 2tt * Re # (45 ) 

c c v ' 


MJkm-tr* MAFtuE-m* co#?poi?yino/v 

DENVER DIVISION 



7-14 


MCR-68-87 


Therefore* the unit radial strain (which is equal to AR divided 
by the undeformed radius) is 

e = AR/R = £ • (46) 

r c 

The unit longitudinal module strain, based upon the equivalent 
longitudinal module elasticity (Section 7.1.2) can be expressed as 

e x = F/AE = F/[2n(d/£) 3 AE)] (47) 

where AE is defined in Equation 8 and AE defines the elasticity of 
an individual truss rod* 


It should be noted that these results are based upon the 
assumption that there is a large number of truss bars so that the 
concentrated radial joint forces can be approximated by a uniformly 
distributed force* The concepts presented here are also applicable 
to other force components such as those which lead to bending and 
to lateral displacements. The extension of the analysis to include 
these effects is omitted here. 


7.1.7 Circumferential Bending Stiffness 

Proper simulation of the shell ring modes requires the simulation 
of circumferential bending stiffness. The bending stiffness of the 
rings can be expressed as 


where 


d 


El = dD 


r r w 

longitudinal distance between ring centerlines* 


(48) 

and 


D = bending stiffness of a circumferential strip of width w 
for the prototype shell. 


7*1*8 Computer Program 

Appendix B presents the listing of a computer program* written 
in FORTRAN IV for the IBM 1130* to compute truss-ring module 
parameters from full-scale shell parameters. 
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7 . 2 Use of Reinforced and Stiffened Shells for Simulation of 

Orthotropic Cylinders 

Although the structure to be modeled is a cylindrical shell) 
the following developments can be carried out with little regard 
to the shape of the structure. Flat plate elements are used and 
all materials involved are assumed to obey the generalized Hook's 
law. Bending and extensions! rigidities for an orthotropic plate 
are found and the equivalent orthotropic properties are expressed 
in terms of the model shell geometry and constituent material 
properties . 


7.2.1 Bending and Extensional Rigidities for an Orthotopic 
Plate 


An orthotropic plate is one which has different material prop- 
erties in two mutually perpendicular directions. In this section 
the two directions will be referred to as the longitudinal and 
transverse. Hooke's law for an orthotropic plate iss 


C 11 e i + C 12 e t 
C 12 e S. + C 22 e t 


(49) 


T Jlt ‘ C 66 Kit 


By applying direct stresses in each of the two directions) the 
following set of relations can be obtained. 
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Equations (49) can be rewritten to give strain instead of stress. 
1 


(C 11 C 22 ” C 12 J 


(C 22 °l~ C 12 V 


(C 11 C 22 “ C 12 J 


<“ C 12 + c u a t ) 


4 T)lt 


or 


^ ■ 




e^ ( * w a + V and Ut ' g V 


It can be seen that there are four independent material properties 
necessary to describe the behavior of an orthotropic plate. These 
properties are G and any three of the properties in Equation (50) . 

The existence of Equation (50) dictates that G must be one of the 
independent properties. 

Consider the bending of an orthotropic plate in the longitudinal 
direction. For small deflections the curvature can be expressed by 



where u is the lateral deflection. The unit elongation a distance 
a from the middle surface is 


- z 


a 2 o) 

a i 2 


In order to maintain continuity during bending the transverse strain 
must be zero. Therefore? 

H ' eJ <0 <i - V a t> 


; t ■ r (-"tA + V - °* ° t ■ “m 
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Therefore} 


and 


ir u - w 


At y tA* 
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The total bending moment can be obtained by integrating the bending 
stresses through the thickness of the plate 


M 


h/2 h/2 

In ’ 'In u ' ^ 3 * 2 " z 
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A d 


E I 
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Therefore, the bending rigidity for longitudinal bending isi 


* 12(1 - VtJ 


or 


2 3 

e t h 


2 12 <*» - ViV 


A similar treatment for bending in the transverse direction yields: 

W s 

D 


1 12 < E * - E t ^n 2 e > 

Expressions for extensional rigidities in each of the two 
directions can be obtained in much the same manner. 
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7,2,2 Equivalent Properties in Terms of the Constituent 
Material Properties and Shell Geometry 
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Flg 9 7.7 Typical Model Section 


Nomenclatures 

u matrix Poisson's ratio 
r m 

E ,G matrix moduli 
m m 

fiber moduli 

fiber cross-sectional area 
h thickness of unreinforced shell 
t depth of circumferential stiffeners 
width of stiffeners 
stiffener spacing 
longitudinal fiber spacing 
transverse (circumferential) fiber spacing 
R middle surface radius 

In order to derive an expression for the equivalent modulus 
in the longitudinal direction, a pseudo material is first developed. 
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This material consists of an unstiffened plate with fibers in the 
transverse direction only. The fibers are assumed to have a square 
cross-section with the same area as the actual fiber. This simpli- 
fying assumption allows one to assume a uniform strain distribution 
on a plane tangent to the fiber and perpendicular to the longitudinal 
axis. Equilibrium is enforced across this plane between adjacent 
elements with and without an encased fiber, see Fig. 7.8. 



Figure 7.8 Unstiffened Plate, Transverse Fibers 


The above conditions and assumptions lead to the following expression 
for the longitudinal modulus for a material with transverse 
fibers only. 


Where 



S t E m UE f + (l-X)Ej 

bE + (S -b)[XE,+(l-X)El 
me i m 



Now, suppose that the stiffeners are included and that their 
presence is accounted for by allowing h to be replaced by an 
effective thickness h, such that equilibrium is maintained between 
stiffened and unstiffened portions of the shell. 


Then: 


and 



h(S x + S 2 )(h + 2t) 

h(S 1 + s”) + 2tS 2 
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where X “ “ 
h 

The total effective longitudinal modulus can now be obtained, 
using as the matrix modulus for longitudinal extension* The 
total load in the 2 . direction will be the sum of the loads carried 
by the longitudinal fibers and the matrix, A uniform strain dis- 
tribution is assumed. 



Let 



then 


- E f a + E^(l - a) 


The derivation for the extensional modulus in the transverse 
direction is similar to the preceeding, however the effective 
thickness, h, is replaced by h*. The new effective thickness is 
obtained by allowing the stiffeners to be "smeared" over the plate 
surface such that an equivalent area is offered to resist trans- 
verse loading. 


h(S 1 + S 2 ) + 2tS 1 

(s x + s 2 ) 


Thens 

S 0 E [A'E + (l-A')E ] 

- x, m f m 

t “ bE + (S p - b)[XE f + (l-A’)E . ] 
m Jc r m 


and 

where 


E t - E f e + E t (l-8) 

^f A 

6 ■ a 7 ■ s^- and x ' 
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It has been demonstrated that the extensional modulus in each 
of the orthogonal directions is a function of the fiber fraction. 
Each material contributes to the equivalent modulus according to 
its percent of volume present. However, an equivalent modulus 
derived from bending considerations depends upon the donation of 
each constituent material according to its percent of inertia 
present. For this reason, two different moduli result, depending 
on the mode of deflection; an extensional modulus and a bending 
modulus. In regard to the rigidities obtained in Section 7.2,1, 
when the product E^I^/I is not a small quantity, the appropriate 
modulus must be used in each rigidity expression according to the 
implied deflection scheme. 


In order to maintain consistency with the derivation of the 
longitudinal extensional modulus, the stiffeners are assumed to 
influence longitudinal bending in accordance with h. Maintenance 
of compatible deflections and slopes of the fibers and matrix 
material during bending yields the following expression for the 
longitudinal bending modulus. Transverse fibers are assumed not 
to influence longitudinal bending. 


hi ■ h y i + 



Transverse bending involves the inertia of the transverse 
fibers and stiffeners. The stiffeners are accounted for by 
using the equivalent thickness, h'. Longitudinal fibers are 
assumed ineffectual during transverse bending. 


ht m E f Y f + E m (1 'V 



Consider the following model, shown in Fig, 7.9. 

For a load applied in the longitudinal direction, there will be a 
corresponding strain in the transverse direction, which is 
proportional to the longitudinal strain. The proportionality 
constant, or Poisson's ratio U^ t , can be obtained by enforcing 
transverse equilibrium between fiber and matrix material and 
assuming a unifrom strain distribution on a plane anywhere in 
the element, provided that it is normal to the transverse axis. 
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Fig. 7 .9 Model of Element Reinforced in Two Directions 


For transverse equilibrium, 

a l A l + a 2 A 2 “ °3 A 3 + °4 A 4 



Fig. 7.10 Notation of an Element 
where « h'S t - b 2 



= (h' - b) S - b 2 
a 4 - b(S fc + b) 

The equilibrium equation becomes? 

e i {E m (h ' S t~ b2 > + E f b2} " e 3 {E m [(h f “b 2 ] + E f [b(S t +b)]} 

With + e 3 " e t 
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or > s ■ 6 t s * ■ e i <s r b) + e 3 b 


and 

Therefore: 


*1 ■ b m e £ 


Hn .. , , , , E „< h ’ S t- b2 > + E f » 2 

Wot " sT < s r b > + b 2 

1 E m [(h'-b)S t -b Z ] + E f [b(S t +b)] 


The method for determining an expression for the equivalent 
shear modulus is analogous to that used for the equivalent 
extensional moduli , A pseudo material is developed for each 
orthogonal direction which includes the effects of resistance to 
shear in the direction of each fiber axis only. A uniform strain 
distribution is assumed on a plane tangent to the fiber and 
perpendicular to the plane, and equilibrium is enforced across 
this plane. For the total equivalent shear modulus the shear moduli 
of the pseudo materials are used as the matrix shear modulus in 
the appropriate direction. See Fig, 7,11 for notation. 



Fig, 7,11 Shears on an Element of Composite 
For equilibrium, 


with 


p t s t - V* 


p t - V* h ’ 


and 


P £ “ W 
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Therefore, 


Also, 

P t - G £ Y t b 2 + Vt <s i h ’- b 2> ■ W 

h - W' + ■ Vt h ' 


The total shear strain is 


Y 


Y*. + Yt 




Therefore, 

2[G f b 2 +G t (S £ h , -b 2 ) ] [G f b 2 *fC A <$ t h'-b 2 ) ] 

S^h’ [G f b 2 +G £ (S t h’-b 2 ) ] + S t h' [G f b 2 +G t (S jl h , -b 2 ) ] 

where, 

S t G m [G f X " +G m (1 ” X " ) 1 

G t “ bG +(S -b) [G,X"+G (1-X")] 
m t r m 

S.G [Gi"+G (W")l 
X, m r m 

G £ " bG m +(S r b)lG f X"+G m (l-X")] 


In retrospect, it should be mentioned that the preceeding 
analysis involves simplifying assumptions which should remain true 
for small displacements of thin plates (or shells). The intent of 
the analysis is to provide relatively simple, but meaningful 
expressions for the equivalent, orthotropic properties. Model 
geometry and constituent material properties are parameterized in 
order to reproduce prescribed equivalent properties. 
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8. APPLICATION TO SPACE STRUCTURES 


The design and development of flexible orbital structures 
proposed for the future will require an accurate prediction of the 
structural deformations and dynamic loads encountered in a space 
environment. Some of the typical configurations that may be 
deployed during the next decade are; 

a. Orbiting space stations such as the Apollo Applications 
Program (AAP) Orbital Workshop and the Manned Orbital 
Laboratory (MOL). The maximum dimension of such stations 
may be several hundred meters; 

b. Dish antennas or radio telescopes with maximum dimensions 
of several hundred meters; very large radio telescopes with 
typical dimensions of several thousand meters. 

The longest characteristic elastic periods of these structures 
may range from approximately 5 seconds for the manned space stations 
to values greater than 1000 seconds for the radio telescopes . The 
dynamic disturbances which influence these large artificial satellites 
are, in many instances, of a magnitude that is insignificant in the 
analysis of current, relatively stiff, space structures. 

The dynamic behavior of AAP and MOL class space stations is 
characterized by the interaction between various dynamic disturbances 
including elastic deformations, rigid body motions and liquid motions. 
The dynamic behavior of large antennas or radio telescope structures, 
deployed and stabilized by centrifugal forces, can be affected by the 
interaction between the geomagnetic field and induced currents. 
Additionally, thermal distortions and distortions induced through 
photon pressure and the interaction of spin rate with the gravity 
gradient may produce appreciable elastic deformations. 

An indicative (but not dimensionless) parameter that may be used 
to define the regimes of dynamic model simulation for vibration testing 
is the logarithm of the quotient of the characteristic maximum prototype 
dimension divided by the square of the fundamental prototype frequency, 
log (l/f2). Table 8.1 shows approximate model simulation regimes. The 
characteristic dimension is measured in meters and the frequency in 
Hertz. 
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Table 8.1 Approximate Model Simulation 
Regimes for Vibration Testing 


Scaling Parameter 
Log a/ f 2 ) 

Method of Dynamic 
Simulation 

Conf iguration 
Applicability 

-2 

0 

2 

Predominately full- 
scale testing; some 
replica scale models 

Small Launch Vehicles 

Direct geometric 

AAP, MOL* etc.; 


models; some distorted 

Large Launch Vehicles 


geometric models 


* 

General geometric 
models (geometric 

100-Meter Antenna* etc. 

6 

8 

10 

simulation* distorted 
physical properties) 


Completely distorted 

1500-Meter Radio Tele- 


models 

scope* etc. 


8 ,1 Scaling of Manned Orbital Configurations 

The use of direct geometric scaling using similar materials is> 
in general, feasible for these configurations. This replica scaling 
will automatically satisfy structural scaling laws but* unfortunately 
is restricted by manufacturing and tolerance control limitations that 
arise in small scale models. Also, replica scaling cannot be used to 
simulate slosh effects. These limitations are discussed in Chapter 9. 

In those areas where direct replica scaling is not possible, equivalent 
distorted scaling techniques must be considered. Distorted scaling of 
shell-like structures is discussed in Chapter 7. The equivalent 
distorted scaling techniques must be carefully analyzed to ensure that 
the portion of the structure simulated by these techniques will deform 
in a manner similar to the prototype. Distorted scaling usually 
neglects some scaling parameters that are not important to the simula- 
tion of an individual phenomenon in order to preserve the more important 
parameters . 
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8.2 Scaling of Antennas and Radio Telescopes 

The overall complexity and relatively small dimensions of the 
individual structural members of these configurations precludes 
the use of replica scaling. Adequate component strength for assembly 
and testing in the 1-g earth environment requires a much stiffer 
structure than would result from direct scaling. Frequencies of 
the model structures typically should exceed by one or two orders 
of magnitude the frequencies which result from direct scaling. 

Model masses may be two to four orders of magnitude higher than the 
equivalent directly scaled masses. In general, it should be possible 
to construct a limited distorted dynamic model that preserves the 
overall geometry but does not reproduce exactly the local geometries 
of the prototype. 

Because the structural frequencies of large orbital structures 
are very low, the elastic vibration modes will couple with the gravity 
gradient induced mode (Ref 8.1). These modes will also respond to 
the excitation caused by local variations of the geocentric gravity 
field. In addition to gravity gradient effects, large structural 
deformations may be caused by periodic thermal expansion and by the 
incident solar wind. 

MacNeal (Ref 8.2) discusses several of these effects in relation 
to a typical 1500-meter radio telescope and concludes that the maximum 
total photon pressure load is approximately 1/15 of the axial component 
of tension in the spiral fibers due to spin. Further, when it is assumed 
that the ratio of fiber stress to axial load is the same for the two 
conditions, the resulting reflector stresses are approximately ten 
times the stress due to the gravity gradient. A comprehensive investi- 
gation of gravity gradient- induced stresses in an axisymmetric body 
rotating about an axis perpendicular to the orbital plane can be 
found in Reference 8.3. 

8.2.1 Simulation of Thermal Effects 


The thermal modeling analysis of space structures is based upon the 
assumption that the structure is operating in a high vacuum environ- 
ment. Therefore, the only thermal considerations are heat transfer by 
radiation and heat conduction within the structure. 

The functional expression for the rate of internal heat production 
is 

q = f (k, H, e, a, £, a, l 9 t„ T) 


M/8/7T//V 

DENVER DIVISION 



8-4 


MCR-68-87 


where 


3 

q = rate of internal heat production/unit volume, Q/L T 

k - thermal conductivity, Q/LTR 

/ 2 

H = radiation intensity, Q/L T 

3 

c = heat capacity/unit volume, Q/L R 

2 4 

O' - Stef an-Boltzmann constant, Q/L TR 
£ = emissivity, 
ot = absorptivity, 

H = characteristic length, L 
t = time, T 
T = temperature, R. 


The basic variable dimensions (Q, L, T, R) yield six 
independent scaling relations 


2 

= e, 7?2 = a, tt 3 = kt/c& , 

tt^ = ax 3 #,/ k, IT = H&/kT, and IT = q& 2 /kT . 

Other sets of non-dimensional groups might be obtained from 
the fundamental expression. This group is particularly convenient 
for dynamic scale modeling as each term, where possible, contains 
the characteristic length parameter (£) . Thus, it is possible to 
relate the several scaling parameters to the mo d el- to- pro to ty p e 
scale factor. 

Two methods of thermal simulation can be developed. These methods 

are; 

a. Maintain identical temperature time histories for model (m) 
and prototype (p) with no material restrictions, or 

b . Maintain identical materials for model (m) and prototype (p) 
and predict the temperature time histories. 
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The similitude requirements for the two cases, developed from 

7T- through TTg , are presented in Table 8,2. The parameter n defines 

the length scale factor, n = 1 / 9 , . 

■ • m p 

Table 8 ,2 Thermal Modeling Similitude Requirements 



The concept of maintaining identical temperature time histories 
appears to be advantageous where steady state temperature conditions 
are involved. The advantage of either method is not obvious when 
transient thermal conditions are considered. A complete discussion of 
thermal similitude requirements for spacecraft models can be found 
in Reference 8.4. 

8.2. 2 Simulation of Gravity Gradient Effects 

The necessity to scale the influence of the gravity gradient 
on large space structures depends upon the oscillatory period of the 
fundamental mode. When this period exceeds about 1000 seconds, the 
orbital structure cannot be regarded as an independent system. At 
higher structural frequencies there is sufficient frequency separa- 
tion between elastic vibrations and gravity gradient-induced rigid- 
body oscillations and the coupling between these two types of modes 
can be disregarded. 
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Consider, for example, a beam-like space structure rotating about 
its center of mass (Fig. 8.1). A circular orbit is assumed. 



Figure 8.1 Beam-Like Space Structure, Circular Orbit 

\ 


The longitudinal structural deformation can be expressed in 
functional form as 

6 = f (m, £, AE, \p 9 3G/3R, t, t) 

where 6 = longitudinal deformation, L 

Z = characteristic length, L 

AE = extensional stiffness, F 

\p = beam rotation rate, 1/T 

3G/3R = gravity gradient, 1/T^ 
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T = orbital period, T 
t = time, T. 

The functional relationship contains eight variables in three 
basic dimensions. The principles of Chapter 3 indicate that five 
non-dimensional scaling relations exist, They are 

= 6/£, tt 2 = m£ 2 ip 2 /AE, = tip, 

tt a = t/T and 7T 5 = (3G/8R)x 2 . 


The scaling relationships obtained from the above can be expressed 
in terms of the length scale (n = #^) and the material properties as 

6 /6 = n, iW^_ = tm (AE) /m (AE) ] 1/,2 /n, 
mp mp p mm p 

t /t = Ip /ip , T /T = ip/'P, 

mp pm mp pm 


and (3G/3R) /(3G/3R) = (\|> /ip ) . 

m p m p 


The use of a centrifuge for model simulation of the gravity 
gradient is considered to be feasible. A possible test configuration 



\ 


/ 


\ 


V 




/ 


Model 


Figure 8.2 Centrifuge for Gravity Gradient Simulation 
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The desired model gravity gradient defines the required 
centrifuge velocity, 

ft = (3G/3R) 1/2 . 

m 

The period scaling, defined by the above expression and is 

T /x = (3G/3R) 1 ^ 2 /0, 
m p p 

and the velocity scaling is 


ip /ty = fl/ (3G/ 3R) V 2 

m p p 

The stiffness scaling, from is 

(AE) /(AE) = n 2 (m /m )ft 2 /(3G/3R) . 

m p m p p 

An additional angular velocity, ft , must be applied initially 
to the model to compensate for the relative motion between model 
and artificial gravity gradient vector. This velocity is 

ft = ft + 1/x . 
m m 
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9. LIMITATIONS AND ECONOMICS OF DYNAMIC MODELS 


The following sections define some of the major limitations 
associated with the simulation of full scale dynamic characteristics 
and the economic aspects of model design and construction* 

9,1 Limitations 


i) Manufacturing limitations— for very small scale factors 
it becomes very difficult to reproduce exactly the local 
geometry of the full scale configuration; errors due to 
material tolerances are increased. 

ii) Non-linear structural effects— at certain loads conditions 
local skin buckling may occur. 

ill) Aerodynamic effects— direct scaling of aerodynamic effects 
may not be possible under certain test conditions. 

9.1,1 Manufacturing Limitations 

The reproduction of launch vehicle configurations to small 
scale factors is restricted by the manufacturing limitations. 

Scale factors as low as 0.30 should cause no manufacturing problems 
for current launch vehicle configurations; at smaller scale it 
becomes increasingly more difficult to reproduce exactly the local 
geometry of stiffened shells, truss-like structures, and propellent 
tank walls. In addition, rivets and other fasteners become 
unreasonably small and the assembly of components by welding is no 
longer possible. An excellent example is found through examination 
of the prototype Saturn V vehicle. Typical second stage tank 
geometry and second-to-third stage adapter geometry (Ref 9,1) are 
shown in Figures 9.1 and 9.2. Note that at 10% scale factor skin 
thicknesses less than 0,005" would be realized in a model that is 
almost 40 * high. 

In order to keep model manufacturing costs within allowable 
limits, it is usually necessary to work with commercially available 
stock and thus eliminate costly manufacturing processes. Tolerances 
on available stock materials limit the accuracy of model calculations 
and the use of such stock introduces increasingly serious limitations 
as the proposed scale factor is diminished. The problem can be 
minimized by selective sizing of stock and by measuring the actual 
dimensions of the scaled structure. 
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Figure 9.1 Tank Wall Properties, Second Stage, Saturn V 



Figure 9,2 Adapter Properties, Second to Third Stage, Saturn V 
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9.1.2 Non-Linear Structural Effects 

Several existing directly scaled models have exhibited a high 
degree of fidelity in their ability to reproduce accurately the 
dynamic characteristics of the prototype vehicle; the two best 
examples being the 10%-scale Saturn V (Ref 9.2) and the 20%-scale 
Titan III (Ref 9.3) models. However , while these models have 
exhibited excellent overall dynamic reproduction of their full- 
scale counterparts, in certain instances they have also exhibited 
certain of the limitations of the direct scaling approach. The 
necessity to reproduce the shell structure which comprises a major 
portion of this type of configuration requires that the model be 
relatively large. Nevertheless, in certain areas of the vehicle 
the resulting skin gages may still be Inadequate to resist the 
assembly stresses caused by the integration of several components 
into the complete model. Imperfections of this type are inevitable 
when direct geometric scaling is employed and it may be necessary 
in the future to consider scaling techniques whereby direct 
scaling is applied to those areas of the vehicle where skin gage 
will permit and other, more critical areas employ distorted 
scaling , 

9.1.3 Aerodynamic Effects 

The direct scaling has only limited application in the area 
of aerodynamic testing due primarily to the restraints imposed by 
available tunnel facilities. The investigation of wind-induced 
oscillations on a Titan III model (Ref 9.4) illustrates this point, 
A 7,5% dynamically scaled model of the Titan III vehicle, trans- 
porter launch frame and umbilical mast was constructed in addition 
to a 7,5% geometrically scaled model of the umbilcal tower. The 
model simulated the mode shape and frequency characteristics of 
the full scale launch vehicle, umbilical mast and transporter. 

The requirements for accurate dynamic scaling of this type 
of model are: 

Matched Reynolds Number; Re * , 

Matched Strouhal Number; S ■ — * , 

and matched model to air density ratio. 
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However , the tunnel imposed conditions (NASA Langley Research 
Center , 16-foot transonic tunnel) which arose due to the use of 
Freon 12 as the test medium were: 

P TUNNEl/ P AIR " 3,7 
^TUNNEL AIR " 0,757 


and it was therefore not possible to scale the stiffness of an 
unfueled vehicle without exceeding the scaled mass distribution. 

A compromise in the form of a mismatch in either Reynolds number 
or mass distribution was necessary and the former was selected as 
being the least severe. It is important to note that the Reynolds 
number of the model test in the tunnel will usually be much lower 
than the true flight condition (in the wind-induced oscillation 
test it was 45% of the full-scale value) but there is little that 
can be done to alleviate this problem. Test conditions should be 
chosen to make the Reynolds number as high as possible but the 
full-scale value will seldom be approached , 

9,1,4 Structural Damping Effects 

Structural damping is a phenomenon strongly dependent not 
only on the type of construction but on material finish, pressure 
in the joints caused by rivets or other fasteners (Ref 9,5), non- 
structural items mounted on primary structure, etc. Limited data 
is available on modal damping of models and for that reason no 
statistical evaluation is possible. However, a rather complete 
set of damping ratios was obtained during the Titan III model 
vibration test (Ref 2.3) . Figure 2,8 (taken from Ref 2,3) shows 
the trend of damping ratios for lateral and longitudinal modes , 

No full scale corresponding values are available, consequently 
no conclusions can be drawn. The inability to assess the prototype 
model damping coefficient from model tests makes the use of 
dynamically similar models for design purposes somewhat question- 
able. This fact becomes even more so when unconventional structures 
are dealt with where no known previous damping data are available, 

9.2 Economic Aspects 

It does not seem logical to design and fabricate a model 
whose overall cost is greater than that of a prototype. On the 
other hand the manufacture of an inexpensive model may result in 
misleading and even incorrect data due to manufacturing simpli- 
fications which yield an incorrect dynamic simulation of the 
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full scale structure* It becomes essential to consider the overall 
optimization of the economic situation in order to achieve a 
realistic model test program - a program yielding results within 
prescribed tolerance levels and, at the same time, a program within 
allowable budgetary limits. 

The study of the economics of dynamic model testing may, in 
a very general sense, be divided into four major categories t 

i) Model design and construction 

ii) Required test facilities, equipment and instrumentation 

iii) Model testing 

iv) Data acquisition and analysis. 

Approximately 30 elastically scaled dynamic models of 
aircraft and launch vehicles were surveyed to establish dynamic 
modeling costs. Data collected included method of construction, 
total cost of model, estimated engineering cost, material cost, 
and a measure of simulation quality (the approximate number of 
modes simulated by the model). It was established that the method 
of construction and the approximate number of modes simulated, 
rather than model size or weight, defined model costs. The summary 
of the results of this survey is shown in Figure 9.3. 
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Figure 9*3 Dynamic Model Cost 
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10 * CONCLUSIONS AND RECOMMENDATIONS 


A comprehensive summary of the dynamic scaling laws applicable 
to current launch vehicle configurations is presented * Primary 
consideration is directed to the study of direct geometric scaling 
using similar materials (replica scaling) although, where applicable* 
certain distorted scaling techniques are considered* 

A review of the current state-of-the-art regarding launch vehicle 
dynamic models indicates that the several vibration survey models 
used to either augment or replace the usual analytical investigation 
have exhibited faithful dynamic representation. Additionally, 
dynamic models used for aeroelastic investigations (e.g., wind- 
induced oscillations, buffet) have yielded satisfactory results in 
view of the limiting conditions imposed by available wind tunnels* 

A more extensive investigation of aerodynamic scaling and the 
restrictions imposed by the failure to exactly match all the required 
scaling parameters is recommended. 

A simple matrix procedure which regularizes the derivation of 
scaling parameters to the extent that they can be developed through 
digital computer analysis is detailed. This procedure has no known 
previous application and is considered to be a substantial improve- 
ment over existing methods* Most of the scaling parameters presented 
here were developed using this procedure. 

The scaling laws applicable to beam-like structures are 
summarized and the advantage of replica scaling, where possible, is 
emphasized. Simulation of uncoupled bending or torsional vibration 
is discussed and it is concluded that replica scaling is, where 
viscoelastic materials are concerned, impossible. An investigation 
of propellant slosh effects in beam-like tanks indicates the 
feasibility of simulating boost acceleration levels in a 1-g test 
environment through the proper adjustment of physical properties 
and time scale. A study of load and response scaling indicates that 
geometric similarity is not required to simulate these phenomena. 

Simulation requirements for shell-like structures including 
isotropic and orthotropic cylindrical shells are summarized* This 
summary includes a study of coupled elastic-slosh modes and it is 
concluded that distorted scaling is required if small scale models 
are used to simulate the boost acceleration environment. The need 
for further investigation of possible construction techniques to 
simulate slosh-elastic coupling is clearly indicated. An experimental 
program to build and test such a model, possibly using plastics 
loaded with metal powder and heavy liquids as the simulated 
propellants is recommended. 
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Scaling laws applicable to the simulation of both uncoupled 
and elastically coupled lateral propellant sloshing are developed 
and it is concluded that an arbitrary prototype acceleration can 
be simulated with an appropriate adjustment of the time scale. 

The difficulty is maintaining exact similitude in all structural 
parameters, due primarily to manufacturing limitations imposed 
upon small scale models, is emphasized. Scaling requirements for 
simultaneous simulation of two or more liquid effects are presented. 

The distorted shell simulation methods developed during this 
study (truss-ring modeling, reinforced plastic modeling) have had 
no previous practical application. An experimental program oriented 
toward verification of these methods is being performed. The 
results of this study will give a reasonable indication of the 
applicability of these distorted scaling techniques. 

Only a limited investigation of modeling techniques applicable 
to large space structures has been conducted. A serious need exists 
for a detailed theoretical and experimental program to investigate 
the modeling techniques required for the simulation of such structures 
in the test laboratory. Simulation of gravity gradient effects, 
photon pressure effects, cyclic radiative thermal effects and 
interaction of induced currents with the Earth’s magnetic field are 
possible areas of study. These effects, while insignificant to the 
"small 11 and rigid structures of the usual Earth environment (including 
launch vehicles), will couple with the structural dynamic deformations 
of large, and relatively light, orbital structures. The modeling 
techniques necessary to simulate these structures in a space 
environment will be established only through continued and intense 
research studies. 
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B.l ORTHOTROPIC SHELL SIMULATION USING TRUSS-RING MODULE 

The scaling parameters for simulation of orthotropic cylindrical 
shells are 


2 2 

d X to /K , 
s xx’ 

*2 

- K /K , 
xy xx’ 

*3 

- K /K 
yy xx 

D /K X 2 , 

XX XX 9 

*5 

- D /K X 2 , 
xy xx 9 

*6 

- D /K 
yy xx 

B /K , 
xy xx s 

*8 

- B /K , 
xz xx’ 

"9 

» B /K 
yz xx 


and “ X/J l 


where w * frequency, 1/T 

2 3 

d g ** surface mass density, FT /L 

K ,K ,L * extensional rigidities, F/L 
xx xy yy 

D ,D ,D * bending rigidities, FL 
xx xy yy 

B » in-plane shear rigidity, F/L 
xy 

B xz ,By z * transverse shear rigidities, F/L 
S, = characteristic length, L 
X ° any length. 

The relations between orthotropic shell parameters and truss 
ring module parameters (Section 7.1,6) are 


Orthotropic Shell 
K 

xx 

K 

yy 


yy 

B 

xy 


Truss-Ring Module 

(n/TT) (AE/R) (d/S,) 3 

A E /d 
r r 

El/d 
r r 

(n/TT) (AE/R) (d/S,) (a/ S,) 2 cos (a/2) 
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K 


yx 


V 

e 


Truss-Ring Module 

(n/tf) (AE/2R) (d / £) (a/ £) 2 

(n/tr) (E/E r ) (A/A r > (a/d) (d/A) 3 (a/2R) 


where n = number of joints on one side of ring, 

R » radius of ring centerline, 

d “ longitudinal distance between ring centerlines, 

2 2 2 

A * length of truss bar, Z * a + d , 
a = ir/n, 

a = 2R sin(a/ 2) , 


A * area of truss bar, 

E = modulus of elasticity of truss bar, 

A y ■ area of ring, 

E^ * modulus of elasticity of ring, 

I ■= area moment of inertia of ring. 

and is the equivalent Poisson's ratio. 

The computational procedure for truss-ring module design is 
based upon knowledge of the prototype rigidities. 


K , K , D and B 
xx ’ yy yy xy 

P P P P 


and specified values of the model-to-prototype ratios 


d /d , X A and to /w . 
s s m p m p 

m p r v 


If values for n, R, E and E are selected, the unknown model 
parameters can be established. The truss rod spacing is 

a * 2R sin (a/2) (B.l-1) 
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where a ■ ir/n. The truss ring spacing can be established from 
the scaling relation 


B /K - tt, « B /K - (a/d) 2 cos (a/2) 
xy xx 7 xy xx v ' ' 

ram P P 

and is 

d - a[cos(a/2)K /B ] 1/2 . (B.l-2) 

xx p x ? p 

The truss rod length is then 

l « [a 2 + d 2 ] 1/2 . (B.l-3) 


The truss rod cross-sectional area can be established from the 
similitude requirement 

d X 2 w 2 /K - tt, - d A 2 oj 2 /K 

s m m xx 1 s p p xx 
m m p r r p 

and the scaling relation 

K « (n/ tt) (AE/R) (d / 1) 3 

Xa 

m 


and is 


A - (R/E) (Tr/n) (£/d) K 


xx 


(B.l-4) 


m 


The truss ring cross-sectional area, established by the similitude 
requirement 

K /K - - K /K 

yy xx 3 yy xx 
p p m m 


and the scaling relation 


K 


yy, 


m 


A E /d, 
r r * 


is 


(d/E ) (K 
r /v yy. 


/K )K 


xx 


XX 


m 


(B.l-5) 


Likewise, the truss ring cross-sectional area, established by the 
similitude requirement 
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yy, 


/K. 


xx 


it. 


D /K 

yy 


xx m 
in 


and the scaling relation 

D » E I /d* 
yy r r * 


is 


I - (d/E ) (X /X ) (D /K )K 
r r m p yy xx xx 

r p m 


The equivalent Poisson ratio 


V - (n/7T) (E/E ) (A/A ) (a/d) (d/ A) J (a/2R) 
e r r 


may be written 


v e - [ (n/TT) (AE/R) (d/Jt) 3 ] [ (d/A r E r ) (a/d) 2 ] /2 


V - (a/d) (K /K )/2 
e xx yy 

ra "m 


(B.l-6) 


V * (a/d) (K /K )/2 
e xx yy 

P /7 P 


(B.l-7) 


This procedure gives exact representation of the similitude 
requirements dictated by ir^, tt. t t and tt^ . Close simulation of 
it 2 will be achieved for realistic structures . 
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♦ IOCS (CARD. TYPE WRITER. KEY BOARD# 11 32 PRINT ER .D ISK 1 
C 

C SCALING PROGRAM FOR SHELL SIMULATION 
C INPUT PARAMETERS 
C 

C SL - MODEL LENGTH / FULL-SCALE LENGTH 
C SF z MODEL FREQUENCY / FULL-SCALE FREQUENCY 

C SM - MODEL SURFACE MASS DENSITY / FULL-SCALE SURFACE MASS DENSITY 

C SXX = KXX FOR FULL-SCALE AT THIS SECTION (SECTION 5.21 

C SYY r KYY FOR FULL-SCALE AT THIS SECTION (SECTION 5.21 

C DYY r DYY FOR FULL-SCALE AT THIS SECTION (SECTION 5.21 

C BXY ~ BXY FOR FULL-SCALE AT THIS SECTION (SECTION 5.21 

C 

C REPEAT INPUT DATA FOR EACH CHANGE OF SECTION PROPERTIES 
C STOP WITH BLANK CARD 
C 

1 READ(2»1001 SL * SF »SM. SXX . SY Y.DY Y* BXY 
100 FORMAT ( 7F 10.01 

C 

1F1SL1 500*500.2 
C 

2 CALL BASK(SL*SF»SM.SXX*SYY*DYY*BXY1 
C 

GO TO I 
C 

500 STOP 
END 


SUBROUTINE BASK ( SL* SF » SM * SX X * SY Y . DY Y • BX Y 1 

INPUT PARAMETERS 

N = NUMBER OF JOINTS (N.GT.2 1 
R r RADIUS OF SHELL 
E = TRUSS MODULUS OF ELASTICITY 
ER = RING MODULUS OF ELASTICITY 

READ( 2* 1001 N*R#E* ER 
100 FORMAT (II 0* 3F 10 .01 
C 

IFIN-215# 10*10 
5 WRITE(3.2001 

2 00 F ORMA T ( 40H NOT ENOUGH JOINTS SPEC IF IEO- ABORT CYCLE! 
GO TO 500 
C 

10 C N - N 

ALF A- 3« 14 15 93 /CN 

A:2.*R*SINt ALFA/2.1 

D -A *SQR T { COS( ALF A /2 . 1 * SX X /B X Y 1 

X L-SQRT ( A ** 2* D* *2 1 

SXX M-SX X * SM *SL* *2*SF**2 

A ROD- ( R *3 .141 59 3* XL ** 3/ (E »CN*0* *3 1 1 *SX XM 
ARING-D*SYY *S XX M/ (ER*SXX1 
X IR IN rO* SL *♦ 2* DY Y* SX XM/ ( ER *SXX I 
SYYM-SXXM*SYY /SXX 
P0IS-0.5* A* *2 *S XX M/ ( S Y YM* □♦ *21 
C 

WRITE (3.2 10 1 
c hdm a ria?u 


■>in 
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WRITE (3.201 1SL.SF .SM. SX X . SY Y »DY Y . BX Y 

201 FORMAK 27H MODEL— TO- PRO TOT YPE RATI0S//6H SL -,F7. 3/6H SF r,F7. 3 
1/6H SM = .F7.3//22H PROTOTYPE RIGIDITIES// 7H KXX =.E14.6/ 

27 H KYY z *£ 14 .6 / 7H OYY = .£1 4.6/7H BX Y = .£14.6/// I 
WRITE 13.202 IN* R *E »ER 

202 FORMATl 38H TRUSS-RING MODULE INPUT PARAMETERS/ /20H NUMBER OF JOI 
1NTS - . 14/ 10H RADIUS = .F6.3/31H TRUSS MODULUS OF ELASTICITY = • 
2E14.6/30H RING MODULUS OF ELASTICITY = .£14.6///) 

WRITE (3 »2 03 I A »D.XL» AROD . AR1 NG »X IR IN 

203 FORMAK 37H TRUSS-RING MODULE OUTPUT P AR AME TERS // SH A =»E14.6/ 

15H D =.E14»6/5H L = »E14.6/1 8H TRUSS ROD AREA = »E14.6/19H TRUSS 
2 RING AREA =*E1 4.6/ 37H TRUSS RING AREA MOMENT OF INERTIA = .Ei 4.6/ 
3//» 

WRITEI3.204 1P0IS 

204 FORMAT ( 28H EQUIVALENT POISSON RATIO =.£14.6/// I 
WRITE (3.205 I 

205 F ORMA T ( 51 H SCALING TERMS SATISFIED = P 1 1 II .PI ( 3 > *P 1 1 6) .P I ( 7 > // 
162H SCALING TERMS NOT SATISFIED = PI (4 1 . PI (5 ) . PI (6 I » PI ( 9 > • PI (1 0> / 
2/39H SCALING TERM APPROX SATISFIED = PI ( 2) ) 

500 CONTINUE 
RETURN 
END 
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B.2 ORTHOTROPIC SHELL SIMULATION USING REINFORCED PLASTIC MODEL 


The computational procedure to determine the reinforced shell 
model properties is based upon the physical model shown in Figure 7.7 
and described in detail in Section 7.2, The notation is that of 
Section 7.2, 

The approximate wall thickness, fiber area and fiber spacing 
are determined in the longitudinal direction (perpendicular to 
the ridges) from the simplified rigidity expressions 

K « E h + A.(E, - E )/S„ (B.2-1) 

2 m f f m 2 

and - y 

12D„ - E a + a/(E, - EJ/S 0 . (B.2-2) 

2 m r r m 2 

1/2 

For assumed square fibers, the fiber width is d ■ and, if 

the fiber spacing is specified as S. - nd, Equations B.2-1 and 
B.2-2 become 

K 9 - Eh + d(E- - E)/n (B.2-3) 

x, m r m 

and o o 

12D. ■ E h J + d J (E, - E )/n. (B.2-4) 

x, m t m 

Similarly, the matrix thickness may be expressed as h=kd and it 
is seen that the constants n and k are restricted to N ^ 1 and 
k > 1; the minimum values (n«l, k=l) specify a solid sheet of the 
fiber material. 

Rewriting Equations B.2-3 and B.2-4 in terms of d, k and n 
yields 

K - d(kE + (E, - E )/n) (B.2-5) 

in r m 

and - o 

12D f - cTCkX + (E, - Ej/n). (B.2-6) 

x, m x m 

Figure B-l shows the physical limitations on d as a function 
h. Parametric lines of k-const and N»const for given values of 
and D^ are indicated. 

If k is set at a constant value and d is eliminated from 
Equations B.2-3 and B.2-6, a cubic equation in n can be obtained. 

This equation is 
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where 




A 2 n 


+ A^fl 


+ A 


A„ 


k \ <l 2 V « 2 - h 3 > 

(E f - \)om ^\ 2 - H h 

36 W E f - V 2 
l 2 V E f - v 3 - 


(B.2-7) 


(B# 2-8) 


Physically realistic solutions are limited to k > 1, If the 
value of n obtained for k*l is less than unity, there exists no 
physical solution for the given values of E and E. and the 
specified values of and D^, Realistic values of n and k may 
be obtained by increasing E., A value of n ^ 1 obtained for 
k > 1 represents a physically possible reinforced shell model. 

Further increasing the value of k will result in a decrease in 
the value of n. The algorithm for the optimization of n and k in 
the initial parameter selection maximizes the product (n-l)(k-l). 

When n and k have been established, Equation B.2-5 can be 
used to compute the fiber width. This parameter is 

d =* nK„/[ (nk-l)E + Ej (B.2-9) 

and it follows that the matrix thickness and longitudinal fiber 
spacing are h = kd and = nd. Note that Equation B.2-9 will 
yield only an approximate value for d as Equations B.2-1 and 
B.2-2 are themselves only approximations. After the approximate 
parameters are determined for the transverse direction, these 
longitudinal parameters can be modified by an iteration process . 

The approximate matrix parameters in the transverse direction 
(direction of the ridges) from the simplified rigidity expressions 

K “ E h + 2E Qt + A,(E, - E )/S_ (B.2-10) 

tm m ■ f f m t 

and - , ? 

12D = (1 - Q)E h + Q(h + 2t) E + A, (E, - E ) /S (B.2-11) 

t m m r r m t 

where 

S = width of "ridges" used for bending reinforcement 
(Fig. 7.7), 

T ■ width of "valleys" and 0 “ S/(S+T), 
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The ridge height is t. If it is assumed that A, is the same in 
both the longitudinal and the transverse directions and if the 
transverse fiber sapcing is written as a multiple of the longitudinal 
fiber spacing, Equations B . 2-10 and B . 2-11 can be rewritten as 

K - E h + 2E Qt + d(E- - E )/(nj) (B.2-12) 

t m m r m 

and ij m « 

12D « (l-Q)E h + Q(h+2t) J E + d <E,-E )/(nj) (B.2-13) 

t m m r m 

with the restrictions that 

nj > 1, 0 ^ Q $ 1 and t > 0, (B.2-14) 

The initial data computation algorithm proceeds by arbitrarily 
setting Q = 0,4. Elimination of the product nj from Equations 
B.2-12 and B.2-13 yields a cubic equation for t 

B 3 t 3 + B 0 t 2 + Bjt + B Q - 0 (B.2-15) 

where 

B- - 8QE 
3 x m 

B ? - 12QE h 

B. » 2QE (3h 2 - d 2 ) 
x m 

B n - d 2 (K - Eh) -12D + Eh 3 , 

U t m t m 

If no positive real root exists, the longitudinal and transverse 
directions should be interchanged. If one positive real root 
exists , 


j - d(E, - E )/(nK - nE h - 2nE Ot) (B.2-17) 

" t m t m m 

If more than one positive real root exists, the resulting values 
of j can be computed from Equation B.2-17, No physical solution 
exists for jn 4 1. 

At this point the computer program presents a choice to the 
operator. The current values of j, jn, t, t/h and Q are printed. 

If more than one set of these variables is available the operator 
may select the most desirable solution or he may elect to change 
the value of Q. If this is done, Equation B.2-15 is recycled and 
new solutions are formulated. Once a "desirable" value of j is 
attained, the transverse fiber spacing (S fc = JS^) is computed and 
the preliminary computation of the shell model parameters is complete. 
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The equations of Section 7.2 are then used to compute the 
actual values of D , K , D and K and the equivalent shear 
modulus, G. The ratios of desired to actual shell rigidities are 
determined for the four controlling rigidities and if these ratios 
are satisfactory, the program is terminated. If the ratios are 
not satisfactory, a operator-controlled iteration process is 
initiated and used until the ratios of desired rigidities to 
actual rigidities are computed to within an allowable unit. 
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♦ IOCS (CARD. TYPE WRITER .KEYBOARD. 11 32 PRI NT ER .DISK 1 
C 1130 PROGRAM TO COMPUTE REINFORCED SHELL PROPERTIES 

DIMENSION XCOF! 4) .COF <4 ).ROOTR( 31. ROOT I <3 1 
1 CONTINUE 
C READ OATA 

READ* 2.9991 EM. EF. XMU. GM . GF 
READ! 2.999) SXX. SYY. DXX • DYY 
999 FORMAT ( 8F 10.01 

WRITEC1.998) EM. EF . XMU. GM. GF 
WRITE <1.998 1 SXX. SYY. DXX. DYY 
998 FORMAT!/. 10 X . 5E15.3) 

8 WRITE <1.991 

99 FORMAT! • DETERMINE H.O. AND SLV) 

C ASSI6N SUBSCRIPT T TO LARGER ONE OF BENDING RIGIDITIES AND THE 
C EXTENSIONAL RIGIDITY ASSOCIATED WITH IT 
IF1DXX-0Y Y) 10. 10.12 
10 DL=DXX 
SL-SX X 
DT=DY¥ 

S T-SY Y 
60 TO 14 
12 DL=DYY 
D T-DX X 
SL-SY Y 
S T-SX X 

C START COMPUTATION LOOP ON INDEX L 
C FIRST MMC 1 =MA N— MACHINE COMMUNICATION EXCHANGE) 

C MMC = WHAT K DESIRED 
14 WRITE <1*1 00) 

100 FORMAT! * TYPE VALUE OF K DESIRED IN FORMAT F10.0M 
R EA 0! 6 . 10 1) XK 
1 01 F0RMATCF10.0) 

C COMPUTE POLYNOMIAL COEFFICIENTS TO GET N 
X COF ! 4 1 = X K* *3 ♦EM* < 1 2. *DL* EM ** 2-SL ** 3) 

X COF ! 31-! EF-EM) *! 36 .*DL*XK* *2*EM* *2-SL* *3 1 
X COF ( 2 1 - ! EF-EM) **2* 36 . *DL *X K*EM 
X COF ( 11-1 2. *DL* < EF-EM)* *3 
X COF ( 3 ) -X COF! 3) /XCOF! 4) 

XCOF! 2>=XC0F! 2) /XCOF! 4) 

X COF < 1 ) -X COF I 1) /XCOF! 4 ) 

X COF< 41-1 . 

C POLYNOMIAL ROOTS 

CALL POLRT! XCOF .COF . 3 . ROO TR .ROO TI »IER) 

IF! IER) 16 .1 8. 16 

C MMC = ERROR PROCEDURE ! IF ANY. MAY QUIT. OR START WITH NEW K) 

16 WRITE < 1 .1 02 )IER 

102 FORMAT! • IER=». 13.* TYPE 1.0 FOR NEW K.2.0 FOR STOP * ) 

READ! 6. 10 1) CHECK 

IF! CHECK- 1.0) 17.14. 17 

17 STOP 

C ROOTS OF THE POLYNOMIAL 
C THESE ARE POSSIBLE VALUES OF N 

18 WRITE (1*1 03) 

1 03 FORMAT! • ROOTS OF THE POLYNOMIAL! N ) ARE*) 

DO 20 1-1.3 

20 WRITE <1.1 04 )ROOTR < I ). ROOT I< I) 

104 F OR MAT (/2X.E1 5. 7. • ♦ ' .El 5. 7. • IV) 

C OPTION TO TRY NEW K . OR PROCEED 

WRITE ( 1 .1 05 ) 

1 05 FORMAT! * IF THIS IS NOT SAT ISFACTORY TYPE 2.. IF YES. 1 .* > 



REA016* 10 1) CHECK 
IF! CHECK- 1. >14.22*14 

C MMC = VALUE OF K PREFERRED. VALUE OF COMPATIBLE REAL N 
22 WRITE(1»106) 

1Q6 FORMAT I • TYPE IN DESIRED K IN F 10. 0*) 

READ(6. 10 1) XK 
WRITE U. 1 08 > 

108 FORMAT! • TYPE IN OESIRED VALUE OF N IN F 10.0 •» 

READ! 6. 10 1> XN 

C COMPUTE AND DISPLAY D.H.S1 

D=! XN*SL> /! (XN* XK-1 .» *EM*EF > 

H~XK*D 
S 1=XN*D 

WRITE < 1 .1 10)D.H*S1 

110 FORMAT!* D r *.E15.7.* H = *.E15.7. * SI - * .E15.7) 
WRITE (1*1 11 > 

111 FORMAT! • THESE ARE THE INITIAL VALUES *> 

C END LOOP ON INDEX L 

C START INDEX T COMPUTATION LOOP 
25 WRITEd.l 121 

112 FORMAT! » COMPUTE T. S2. AND 0*1 
C MMC - WHAT IS DESIRED Q 

WRITEil. 11 3) 

113 FORMAT! • SPECIFY INITIAL VALUE OF Q BETWEEN 0 AND 1 .* 
READ! 6. 10 1)0 

C COMPUTE POLYNOMIAL COEFFICIENTS 
XC0F(4I=8.*Q*EM 
X COF ( 31 = 12.*0*EM*H 
X COF ( 2 > =6 *0 *EM* H* *2 -2 * * D* *2 *G *E M 
X COF ( 1 I -D ** 2* ( ST— EM *H 1-12 DT ♦ EM*H**3 

X COF ( 1 ) -X COF ( 1) /X COE 1 4 1 
XCOF! 2) =X COF! 2> /XCOF! 4) 

X COF ( 3>-X COF! 31 /XC0F(4) 

X COF ! 41-1 . 

C ERROR PROCEDURE 

CALL POLRT! XCOF .COF . 3 » ROOTR .ROOTI *IER> 

C IF IER N.E. 0. TRY NEW 0 
C MMC- MAY QUIT ON ERROR. OR TRY NEW Q 
IF! IER >28 .30.28 

28 WRITE (1.1 15 TIER 

115 FORMAT!* IER- * ♦ 13 » * TYPE 1- FOR NEW Q. 2. FOR STOP* 1 
READ! 6. 10 II CHECK 
IF(CHECK— 1«>29»25.29 

29 STOP 

30 CONTINUE 

C POLYNOMIAL ROOTS 
C POSSIBLE VALUES OF T 
WRITE (1.203) 

203 FORMAT! • ROOTS OF THE POLYNOMIAL (T> ARE* I 
DO 31 1-1.3 

31 WRITE (1 .104 IROOTRCI I. ROOTI! II 

C MMC = OPTION TO TRY NEW 0. OR PROCEED 
WRITE (1.1 05 ) 

READ! 6. 10 1) CHECK 
IF ( CHECK — 1« 125.32.25 

C MMC - ENTER COMPATIBLE PAIR OF T AND Q 

32 WRITE (1 .1 171 

117 FORMAT! • TYPE IN VALUE OF 0 * > 

READ! 6* 10 1) Q 

WRITE ( 1 .1 18 > 

1 18 FORMAT! * TYPE IN VALUE OF T *> 
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READS 6t 10 11 T 

X JrD* SEF-EM 1/ S XN*ST-XN*EM *H-2 .* XN *Q *E M* T I 
C DISPLAY VALUE OF J AND NJ 

WRITE II *1 19 1X4 

1 19 FORMATS* THE VALUE OF J IS* »E 14 .7 1 
SAMrX J*XN 
WRITES!. 1201SAM 

1 20 FORMATS • THE PRODUCT J*N IS * .F 10-4 1 
S 2-SI *XJ 

C MC CHOICE TO START T LOOP AGAIN* OR PROCEEO 
WRITE (1 .1 21 1 

1 21 FORMATS * IF THIS IS OK. TYPE 1. « IF NOT 2.* 1 
READS6. 10 11 CHECK 

I F S CHECK— 1. 125* 35.25 

C DISPLAY SUMMARY OF DATA OBTAINEO THUS FAR 
700 WRITES1.1 501 

150 FORMATS • FOR RERUN SPECIFY OESIREO H. D* SI »S2. 0. T IN F10.0M 
READS6. 10 11H 
REAOS 6. 10 1) 0 
READS 6* 10 1) SI 
READS 6. 10 11 S2 
R EADS 6. 10 11 Q 
READS 6. 101) T 

35 WRITES1 .1 24 I H .0 » S 1* S2 . Q » T 

124 FORMATS* SUMMARY OF MATRIX DATA*/* H = *»E14.7/* D z*.E14*7/* Slr**E 
114.7.* S2 -* »E 14 • 7 / • 0 =*.E14.7.* T r*.E14.7> 

C SHOW RATIOS OF DESIRED AND ACTUAL RIGIDITIES 
WRITES! * 1 25 I 

125 FORMATS* SUMMARY OF BASIC RIGIDITY PARAMETERS AND RATIOS TO DESIRE 
ID VALUES* /) 

HBARrl./S H+2. *T 1+2. *T*S 1. -Q )/ IH*S H*2.*T) ) 

XLAMrD*HB AR 

ELB = SS2*EM* SXLAM*EF + S 1.-XLAM1 *E Ml 1 / S D*EM+ < S2-D) * I XL AM*EF* ( 1 .-XL AM ) 
1* EM ) 1 

ALF -D **2*HBAR/S1 
EL-EF * ALF *EL8*S l.-ALF 1 
HPR-1./SH*2.*T*Q1 
X LP -D *HPR 

ET B-Sl *EM* SXLP*EF* SI .-XLP1 *EM1 / S 0* EM* S SI -D 1 * S XLP *EF* ( 1 .-XLP > *EM) 1 
BET -D ** 2* HPR/ S2 
E T -EF *BET *ETB*( l.-BET 1 

XMULT-SXMU/Sl 1* SS1-D)*D*SEM*SS2 /HPR-D** 2 1 *EF* 0* *2 1/ SEM* SSI. /HPR-D 1 
1*S2 -0**21 ♦£ F* ( D * S S2 *0 1 1 1 
C SL -EL * *2 *H/S EL-ET* XHULT* *2 ) 

CSL=CSL/I H* HBAR 1 

CST=EL*ET *H/S EL-ET*XMULT* *2 1 

C ST -CST / ( H* HPR 1 

GAML-3.*D**4*HBAR**3/SS1*3. 14 16 1 
GAMT-3.*D**4*HPR* *3 / S S2 *3 .1 41 6) 

E BL -E F *G A ML *EM* SI .-GAME 1 
E BT -E F*GA MT +EM* S 1 6AMT 1 

CDL -EBL** 2*H* *3/1 12 .* S EBL-E 8T *X MULT ** 2) I 
C DL -CDL / S H* *3 *HBAR* *3 1 

CDT-EBL*EBT *H** 3/ SI 2. * SEBL-EBT* XMULT* *2 1 1 
CDT -CDT/1H**3*HPR**3) 

AF-D* *2 

C TEMPORARY CHECK SECTION 
C ALL DATSWS 3. JJ I 
IF C JJ-21670.6S7.670 
6 67 WRI TE ( 1 .1 75 1 

1 75 FORMATS • SIMPLIFIED PROPERTY CALCULATION* 1 
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CSL-AF*(E F— EM I / SI ♦H*!H*2. *T l*EM/l (H*2.*T) *1 l.-Q »*H*G> 

C ST - AF * !EF-EM)/S2*(H*2*3*T> *EM 

C 

1DL-D* *»♦*! EF— EM) /I 12 .*S1 ) + H* *3*1 H*2. *7 )* *3*EM/ < 1 2. *Q *H** 3* 12 .* < 1 .-0 
l> *!H*2.*T ) * *3 » 

C DT rD ***** (EF— EM 1/ ( 1 2*S2J + 0*EM*I H* 2* T ) ** 3/ 12 ♦( 1-Q 1 *H ** 3*EM/1 2 
670 CONTINUE 

C TEMPORARY CHECK SECTION ENDS HERE 

X SL-CSL/SL 
X ST -CST /S T 
XDL=CDL/DL 
X DT -CDT /DT 

WRITE (1*126 ICSL *SL* XSL 
WRITE (1 .1 27 ICST *ST# XST 
WRITE (1 *1 28 ICDL.DL. XDL 
WRITE(1.129»CDT*DT. XOT 

126 FORMAT!* CSL = **E14.7» * SL -* .El‘t.7.’ RATIO r*.F10.4| 

127 FORMAT! • CST z*,E14 ,7.* ST -• *E 14 .7 #• RATIO = *.F10.4) 

128 FORMAT!* CDL -* .E14 .7** DL -* *E 14 .7 » * RATIO -•»F1D.4I 

1 29 FORMATC* CDT =».».£ 14 *7.. • DT -» *E 14 .7 .* RATIO -'.FI 0.4) 

C MMC r MAY PROCEED. OR RETURN TO T LOOP* OR RETURN TO START 
WRITEd.I 301 

130 FORMA T ( • IF OK. TYPE 0. . IF RESTART • TYPE l.t IF RESTART G* 2.* ) 

R EAD ( 6 » 10 1) CHECK 

IF( CHECK-1. ISO. 8. 25 
C COMPUTE SHEAR RIGIDITY 

50 6-2.* IGF*D* *2 *GM* ( S 1/HPR-D* *2 ) ) 

G=G*!GF*D ** 2* GM * ( S2 /HPR— D ** 2) ) 

6=6/! SI /HPR *( GF *D ** 2+GM* C S2 /HPR-D ** 2) > +S2/HPR *! GF *0 ** 2+GM * ! SI /HPR— 
ID** 2) ) > 

C SS -6 /HPR 
CSLT-CST*XMULT 
CSTL-CSLT *HPR /HBAR 
CDTL-CDT* XMULT 
CDLTrCDTL*! HPR/ HBAR)* *3 
WRITE (1 .1 31 ) 

WRITEd.I 32 ) 

WRITE (1 .1 33 ICSS.CSLT.CSTL »CDTL» COLT 

1 31 FORMAT! • SUMMARY OF THE OTHER RIGIDITIES* ) 

1 32 FORMAT! • CSS CSLT CSTL COTL 

1 COLT*/} 

1 33 FORMAT i 5E 17.8 ) 

WRITE (1 .1 30 ) 

READ! 6. 10 1) CHECK 
IFtCHECK-1. >500.8.25 

500 CONTINUE 
WRITE (1 .301 ) 

301 FORMAT! • TYPE 0 FOR ARBITRARY INPUT. 1. FOR RESTART. 2 TO STOP*) 
READ! 6. 10 1) CHECK 
IF ( CHECK — 1. 1700.1.501 

501 STOP 


SUBROUTINE POLR T ! XCOF »COF »M .ROO TR .ROOT I . IER ) 
DIMENSION XCOF! 1) .C0F<1 l.ROOTR! 1) .ROOT I !1 J 
I F I T-Q 
N -M 
I ER -0 

I F ( XCOF I N *1 ) ) 1 0. 25 .10 

10 I F ( N ) 15.15.32 
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C SET ERROR CODE TO 1 

15 IER-1 
20 RETURN 

C SET ERROR CODE TO 4 

25 1ER-4 

60 TO 20 

C SET ERROR CODE TO 2 

30 I ER-2 
60 TO 20 

32 1 FI N— 36 ) 35.35.30 
35 NX = N 

NXX=N*1 

N2=l 

K J1 = N + l 
00 40 L=I »KJ1 
MT=K J1-L-* 1 
40 COFCMT)=XCOFCL> 

C SET INITIAL VALUES 

45 X 0=1*00500 10 1 

Y 0=0*0100 02 01 

C ZERO INITIAL VALUE COUNTER 

IN-0 
50 X -XO 

C INCREMENT INITIAL VALUES AND COUNTER 

X0--10.0* YO 
Y0=-10.0*X 

C SET X AND Y TO CURRENT VALUE 

X=XO 

Y =Y0 
IN=IN«-1 
60 TO 59 

55 IFITzl 
XPR=X 

Y PR -Y 

C EVALUATE POLYNOMIAL AND DERIVATIVES 

59 I CT -0 

60 UX = 0» 0 
U Y- 0. 0 

Y =0.0 
YT=0.0 
X T= 1 • 0 

U =C OF C N ♦ 1 ) 

IFCU» 65. 130.65 
65 DO 70 1=1 .N 
L = N- 1 + 1 
XT2 = X*XT~Y* YT 
YT2=X*YT.Y*XT 
U =U «-COF CL )*XT2 

Y =V .COF CL l*YT2 
F 1= I 

U X=UX «-F I* XT *COF CL ) 

U Y-UY—F l * YT »C OF CL ) 

X T=XT2 
70 Y T= YT 2 

S UMSO -UX* UX *UY*UY 
1FCSUMS0) 75.110.75 
75 D X- C V *U Y-U* UX )/SUMSO 
X = X«-DX 

DY=-CU*UY *V*UX) /SUM SO 

Y =Y ♦DY 

78 IFC ABS(DY) ♦ ABS I DX »- 1 . OE -0 5 > 1 00.80*80 
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C STEP ITERATION COUNTER 


80 

ICT 

-ICT *1 





I F ( 

ICT 

-50 01 6 0. 

85 

*85 

85 

IFI 

IFI 

T» 100i 

if 90 

.1 

00 

90 

IFC 

IN- 

5) 50.9 5. 

95 




SET 

ERROR 

COOE 

TO 

95 

1ER 

-3 






GO 

TO 

20 




100 

00 

105 

L- 1* NX X 




MT- 

KJ1 

— L* 1 





TEMPrXCOF CMT) 

XCOFCHT I - COFC L 1 
1 05 COF (LI-TEMP 
I TEMP-N 
N -N X 

NXrlTEHP 

IF ( IF I T ) 120*55.120 

110 IF ( IF IT > 115*50.115 

1 15 X-XPR 
Y -Y PR 

120 IFIT-0 

1 22 I F ( ABSIY/X l-l.QE-04) 135*125*1 25 
125 ALPHA-X *X 

SUMSQ-X*X*-Y*Y 
N -N— 2 
GO TO 140 
130 XrO.O 

NXZNX-1 
NXXrNXX-1 
135 Y -0 *0 

SUMSQ-0.0 

ALPHA-X 

N-N-l 

140 C OF C2 > -COFC 2> ♦ALPHA *C0F(1 ) 

145 DO 150 L- 2* N 

150 COFCL^l )-COFf L+ll ♦ALPHA*COF CL l-SUMSQ*COF i L- 1 1 
155 ROOTI (N2) -Y 
ROOTR (N2> -X 
N 2-N2 ♦ 1 

IFCSUMSQ) 160*185*160 
160 Y -— Y 

SUMSQ-0 .0 
GO TO 155 

165 IFtNl 20* 20 .4 5 
END 
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Table C-2 Properties of Aluminum 


FABRI CAB ILIT Y 

WELDABILITY 

HEAT TREATABILITY 

1100, 3003, 5456, 6061 Good 
2014 Fair 
2024 Poor 
7075, 7178 No 

1100, 3003, 5456 Not Hardenable 
2014, 2024, 6061 

Hardenable 

7075, 7178 


ALLOY 

F , KSI 
tu’ 

1100, 3003 

15-20 

2014, 2024 

67 

5456 

50 

6061 

45 

7075, 7178 

75 


FORM 

SIZE, INCHES UNLESS SPECIFIED 

FORGINGS 

MINIMUM MAXIMUM 

THICKNESS 0.100 12 .0 

WIDTH 12.0 +0.063, -0.125 

LENGTH As specified +0.125 50.0 ±1.0 

ROLLED, DRAWN 
BAR 

ROD BAR 

MIN MAX MIN MAX 

THICKNESS 0.375 ± 0.0015 8.0 ±0.031 0.375 ± 0.002 4.0 ± 0,020 

WIDTH 0.375 ± 0,002 10.0 ± 0.062 

LENGTH 12 .0 ± 0,125 50.0 ± 1.0 12.0 ± 0.125 50,0 ± 1.0 

SHEET , PLATE 

SHEET PLATE 

MIN MAX MIN MAX 

THICKNESS 0.006 ± 0.0015 0.249 ± 0,01.1 0.250 ± 0.0025 6.0 ± 0,180 

WIDTH 12.0 ± 0.062 168.0 ± 0.25 12.0 ± 0,375 168.0 ±0.75 

LENGTH 600 .0 ± 1.0 600 ± 1,0 

WIRE 

MIN MAX 

THICKNESS 0.010 ±0.0005 0.375 ± 0.0015 

LENGTH As Specified ± 0.125 5000 ft ± 1.0 

EXTRUSIONS^ 

MIN MAX 

THICKNESS 0 .050 ± 0 .006 7 .0 ± 0 .081 

RADII 0.125 

AREA 0 .2 50 60 .0 

WIDTH 30 .0 ± 0 .081 

LENGTH 500 ft ± 100 ft 

TUBING, PIPE 

MIN MAX 

THICKNESS (WALL) 0.010 ± 0.002 0.500 ± 0.020 

DIAMETER (0 .D .) 0.125 ± 0.003 12,0 ± 0.025 

LENGTH 500.0 ± 100 .0 

^^Mus.t be grain size controlled for better weldability. 
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Table C-3 Properties of Titanium 


FORM 

ALLOYS 

SIZES, INCHES UNLESS SPECIFIED 

FABRICABILITY 

FORGINGS 

COMM. PURE 

CONVENTIONAL FORGING 


MACHINABLE^ 



Ti-6Ai~4V 


LARGE 

SMALL 



Ti-5A^-2 ,5Sn 

THICKNESS 

TOL , +0.25 -0.03 +0 

.09 -0.03 

REQUIRES SCALE REMOVAL 


OTHERS 

LENGTH TOL 

±0.03 

±0.03 

WELDABLE 




MIN . FILLET 


FORMABLE 




RADIUS 

1.0 

0.75 

CHEM. MILLABLE 




MIN. WEB THICK 0.62 

0.25 

HEAT TREATABLE 


ROLLED OR 

ALL^ 


MIN 

MAX 

MACHINABLE 


DRAWN BAR 


1 THICKNESS 

(ROUNDS) 0.3125 + 0.005 1 

WELDABLE 




THICKNESS 

(SQUARE, 0.3125 

4.50 

FORMABLE 




HEX) 



HEAT TREATABLE 




THICKNESS 

(RECT) 0 .200 

10.0 

CHEM. MILLABLE 




LENGTH 


50 ft 



SHEET 

AIL^ 1 ^ 


MIN 

MAX 

SAME AS ROLLED 

OR 



THICKNESS 

0.008 ± 0.002 0.1875 

» ± 0.014 

DRAWN BAR 




WIDTH 

24.0 ± 0.062 100.0 

± 0.125 





LENGTH 

600 .0 

± 0 . 50 



PLATE 

ALL^ 


MIN MAX 

SAME AS ROLLED 

OR 



THICKNESS 

0.1875 ± 0.046 3.0 

±0.10 

DRAWN BAR 




WIDTH 

10.0 ± 0.125 150.0 ± 0.25 






600.0 + 0 .25 



WIRE 

ALL^ 1 ^ 

THICKNESS 

0.009 (MIN. DIAMETER) 

SAME AS ROLLED 

OR 



LENGTH 

500 ft 


DRAWN BAR 


EXTRUSIONS 

COMM. PURE 


MIN 

MAX 

SAME AS ROLLED 

OR 

i 

i 

Ti-6Ai-4V 

THICKNESS 

0 .125 

1 .25 

DRAWN BAR, BUT 

MUST 


Ti-6A,0-6V-2Sn 

FINISH MACHINED 0.060 ± 0.005 


BE DESCALED 



Ti-5Ai-2 . 5Sn 

LENGTH 


40 ft 




Ti-8A,0-lMo-lV 

FILLET RADII 0.062 





Ti - 1 3 V- 1 1 Cr -3 Ai 

LENGTH/THICKNESS 

14 



TUBING, PIPE 

I 

COMM. PURE 


MIN 

MAX 




Ti-6Ai-4V 

WALL THICKNESS 0.004 ± 0.004 0. 

250 ± 0.0250 

SAME AS ROLLED 

OR 


Ti~5Ai~2 .5Sn 
Ti-13V-llCr-3A^ 

DIAMETER(O.D.) 0.250 ^ 12 -° ’o30 

DRAWN BAR 



Ti~8Ai**lMo~lV 

Ti-5Ai-5Sn-5Zr 







^ Unalloyed (Comm. Pure): Ti-0 ,2Pd, Ti-5A.g-2.5Sn, Ti-7A.g-2Cb-lTa , Ti-8Ai-lMo-lV, Tt-6Ai-4V, and 

Ti-6A.g-6V-2Sn. 

(^Approximately 20 (AISI B1112 Steel = 100) . 
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Table C-4 Properties of Nickel 


FORM 

ALLOYS 

SIZES, INCHES UNLESS 

SPECIFIED 

FORGINGS 

INCONEL 600, 700, 


MIN 

MAX 


718 x 750, 901; 

THICKNESS 

0.125 ± 0.005 

4.00 


UDIMET 500, 700; 

DRAFT ANGLE 

2.0 deg 

3 .0 deg 


OTHERS 




ROLLED, OR 

INCONEL 600-702, 


MIN 

MAX 

DRAWN BAR 

718 x 750, 901; 

THICKNESS (ROUNDS) 0.312 - 0.002 

3.00 - 0.005 


HASTELLOY C, R-235; 

THICKNESS (HEX, 




OTHERS 

SQUARE, OCT.) 

0.312 - 0.002 

3.00 - 0.007 



LENGTH 

As Specified 

Spec. +0.025 




±0.125 


SHEET 

ALL ABOVE PLUS 


MIN 

MAX 


WASPALLOY, RENE 41 

THICKNESS 

0.008 ± 0.002 

0.1875 ± 0.014 


NEMONIC 80A, 90; 

WIDTH 

24.0 ± 0.062 

100.0 ± 0.25 


OTHERS 

LENGTH 


600 .0 ± 0 .50 

PLATE 

ALL ABOVE PLUS 


MIN 

MAX 


WASPALLOY, RENE 41 

THICKNESS 

0.1875 + 0.046 

3.0 ± 0.1 


NEMONIC 80A, 90; 

WIDTH 

10.0 ± 0.25 

150.0 ± 0.25 


OTHERS 

LENGTH 


600.0 ± 0 .25 

WIRE 

ALL ABOVE 


MIN 

MAX 



DIAMETER 

0.062 ± 0.0007 

0.500 ± 0.0015 



LENGTH 

As Specified 

As Specified 

EXTRUSIONS 

ALL ABOVE 


MIN 

MAX 



CIRCUMSCRIBED 





CIRCLE 


4.50 



LENGTH 


60 ft 



CROSS SECTION 

0.50 in. 2 




THICKNESS 

0.156 




RADII, CORNER 

0.062 




RADII, FILLET 

0 .250 ± 0 .062 


TUBING, PIPE 

ALL ABOVE 


MIN 

MAX 



WALL THICKNESS 

0 .0015 ± 0.00015 

0.625 ± 0.0625 



DIAMETER (O.D.) 

0.012 ± 0.002 

6.00 ±0 .005 



LENGTH 

As Specified 

100 ft ± 0.25 


Afxm-rgiv MAn imrrA co/?/>o/r>«r/ow 

DENVER DIVISION 







